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SUMMARY 

The  need  to  examine  the  very  foundations  of  fracture  mechanics  became  evident  at 
the  close  of  an  earlier  AFOSR  program  [  1  ].  In  essence,  the  basic  tenet  of  fracture 
mechanics  is  that  the  stress  intensity  factor,  K,  controls  fracture  and  fatigue:  here  the 
former  claim  is  critically  examined.  The  underlying  supporting  arguments  -  the  original 
energy  argument  of  Griffith  and  the  more  rhodern  K-controlled  region  view  -  are 
considered.  These  considerations  demonstrate  that  there  are  questionable  assumptions  in 
both,  so  that  the  viability  of  K  as  a  damage  parameter  for  fracture  has  to  be  established 
by  the  physical  evidence.  The  first  question  then  is  whether  or  not  the  critical  value  of  K, 
K|c>  is  a  material  parameter,  checking  data  shows  it  need  not  be.  The  second  question  is 
can  the  technology  be  usefully  predictive,  even  in  the  most  simple  of  situations:  checking 
the  data  shows  it  to  be  unreliable  in  this  role.  At  this  time  then,  it  remains  to  ask  similar 
questions  concerning  the  role  of  K  in  fatigue  crack  growth  and,  certainly  for  the  monotonic 
loading  case,  develop  alternatives:  these  are  the  objectives  of  the  second  phase  of  the 


program 


INTRODUCTION 

In  general,  the  primary  objective  of  fracture  mechanics  is  to  predict  when  a 
component  will  fail  when  the  analysis  of  the  associated  configuration  leads  to  singular 
stresses.  This  is  not  an  easy  task,  since  the  theory  is  saying  that  physically  unsustainable 
infinite  stresses  result  from  even  infinitesimal  loadings.  At  this  time,  the  accepted  approach 
for  dealing  with  such  anomalous  findings  is  Linear  Elastic  Fracture  Mechanics  (LEFM).  at 
least  in  the  instance  of  cracked  geometries.  Essentially,  LEFM  selects  the  coefficient  of 
the  stress  singularity,  the  stress  intensity  factor,  as  the  parameter  which  governs  fracture 
provided  response  is  sufficiently  brittle.  As  a  result  of  an  earlier  AFOSR  research  program 
[1],  it  became  quite  clear  that  this  choice  and  the  surrounding  technology  needed  a  full 
examination  This  report  describes  work  undertaken  in  a  one  year  study  to  this  end. 

In  somewhat  greater  detail,  the  principal  current  activities  within  a  fracture  mechanics 
treatment  of  a  monotonic  loading  situation  may  be  classified  as  follows: 

•  Singularity  recognition 

•  Evaluation  of  singularity  participation 

•  Interpretation  of  the  resulting  quantified  singular  behaviour 

We  discuss  findings  with  respect  to  each  of  these  in  turn  in  what  follows. 


OVERVIEW  OF  RESULTS  OBTAINED  IN  RESEARCH  PROGRAM 


Singularity  recognition 

Typically  the  presence  and  character  of  singularities  is  well  understood  -  seldom  does 
the  fracture  mechanics  community  stumble  on  this  all-important  first  step.  This  is 
especially  so  for  the  elastic  analysis  of  geometries  containing  cracks.  It  is  also  usually  the 
case  for  the  elastic  treatment  of  other  configurations,  and  for  the  handling  of  cracks  within 
nonlinear  theories.  Asymptotic  analyses  continue  to  be  contributed  in  these  latter  instances, 
eg,  respectively,  [2]"  and  some  excellent  articles  in  the  recent  conference  [3].  These 
last,  incidentally,  demonstrate  that  the  nonphysical  nature  of  singular  stresses  is  retained 
when  one  entertains  yielding,  Le..  the  problem  of  interpreting  singularities  is  not 
circumvented  by  admitting  plastic  flow.  Basically  this  is  because,  physically,  elastic  response 
must  precede  plastic  and  elastic  analyses  to  date  continue  to  be  plagued  by  singularities** 

Singularity  participation 

in  two  dimensions,  analytical  solutions  provide  explicit  and  exact  values  of  singularity 
participation  or  stress  intensity,  K,  in  some  situations.  In  practice,  though,  even  two- 
dimensional  geometries  are  often  too  complex  to  be  tractable  to  a  purely  analytical 
approach.  Here,  however,  numerical  analysis  can  be  employed  to  good  effect  Principally, 
the  finite  element  method,  when  used  in  conjunction  with  a  path  independent  integral, 
furnishes  a  reliable  and  accurate  means  of  determining  K.  This  is  primarily  because  the  finite 
element  method  readily  adapts  to  the  various  and  diverse  geometries  encountered  in 
engineering,  while  path  independent  integrals  are  orthogonal  to  any  other  regular  fields 
present  in  the  analysis  of  a  particular  configuration.  Moreover,  the  approach  is  sufficiently 
efficient  so  as  to  meet  engineering  accuracy  requirements  in  return  for  modest 
computational  effort  Such  efficiency  levels  can  be  demonstrated  on  test  problems  with 
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There  do  exist  elastic  analyses  of  cracks  which  are  apparently  free  from  singularities  -  these,  however, 
either  fail  to  load  the  crack  tip  or  involve  an  integration/finite  length  scale  in  effect  end  therefore  do  not  really 
overcome  the  difficulty. 


known  exact  solutions.  For  example,  the  path  independent  integrals  devised  in  [4]*,  [5]* 

can  be  used  to  determine  K  to  within  1%  using  fewer  than  200  constant-strain-triangle 
elements.  In  all  then,  the  development  of  the  ability  to  establish  singularity  participation  is 
well  in  hand,  at  least  in  two-dimensional  instances.  And  the  extension  of  such  methods  to 
three  dimensions  is  not  seen  as  a  further  development  of  any  great  significance  in  the 
overall  question  of  the  performance  of  fracture  mechanics,  [6]** 

There  are  in  use  at  present,  nonetheless,  some  other  methods  of  determining 

singularity  participation  which  can  be  shown  to  be  quite  unreliable,  [7]*,  [8]*.  In 
substance,  these  methods  employ  some  kind  of  fit  of  local  field  data  to  infer  K:  the  basic 
reasons  for  the  potential  unreliability  of  such  local  procedures  are  as  follows.  First,  all 
local  procedures  must  consider  quantities  near  but  not  at  the  crack  tip.  Second,  at  such 
stations  fields  other  than  the  key  singular  ones  can  contribute  Third,  the  extent  of  such 
participation  cannot  be  either  completely  controlled  or  fully  accounted  for.  As  a 
consequence,  for  any  given  local  procedure  there  exist  problems  on  which  it  produces 
unacceptably  erroneous  results 

A  corollary  of  the  unreliability  of  local  fitting  procedures  is  that  fracture  criteria 
based  on  some  local  field  value  cannot  be  reliabty  connected  to  the  stress  intensity 
involved.  This  means  that  criteria  entailing  crack  opening  displacement  crack  opening  angle, 
or  some  so-called  critical  stress  or  strain  at  a  nearby  station,  all  fail  in  the  second  aspect 
of  fracture  mechanics,  namely  establishing  singularity  participation.  The  successful 
completion  of  this  aspect  is  a  prerequisite  to  successful  interpretatioa  Accordingly  it  is 
unreasonable  to  expect  that  any  of  today  s  local  fracture  criteria  should  prove  of  genuine 
and  significant  value  in  practice 
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**Copy  *ent  earlier  to  AFOSR  (final  report  on  Contract  No  AFOSR-79-007B) 


Singularity  interpratation 

The  first  argument  suggested  in  support  of  the  singularity  coefficient,  or  stress 
intensity  factor,  as  the  controlling  parameter  in  brittle  fracture  was,  in  effect,  Griffiths 
energy  balance.  In  Griffith's  approach,  the  energy  source  is  the  strain  (potential)  energy 
released  on  fracture,  while  the  sink  is  a  surface  energy  term.  Thus  the  controlling  function 
at  the  onset  of  fracture  is  the  energy  release  rate,  G.  Irwin  showed  that  G  has  a  simple 
relationship  to  K  independent  of  geometry,  so  that  Griffith's  energy  balance  is  equivalent  to 
the  choice  of  the  stress  intensity  factor  as  a  damage  parameter. 

Griffith's  work  represented  a  major  step  forward  at  that  time.  Since,  however,  a 
number  of  questions  have  been  raised  concerning  the  extent  of  its  validity  (eg,  Goodier 
[9]).  Moreover,  if  one  simply  applies  the  ideas  to  the  fracture  of  a  truly  brittle  material  in 
a  unaxiat  tension  test,  one  obtains  [10]* 

1 

°°vz  •  (1) 

Here  a u  is  the  ultimate  stress,  l  is  the  length  of  the  specimen.  Equation  (1)  stems  from 
the  fact  that  the  energy  source  has  dimensions  of  F.L,  where  as  the  energy  sink  has 
dimensions  of  F  alone,  F  being  force,  L  being  length.  For  the  application  to  the  tensile 
test.  L  transpires  to  be  t,  the  specimen  length  While  (1)  is  trendwise  correct  in  that 
fracture  stress  decreases  with  increasing  size,  it  is  oversimplified.  This  is  because  a ^  can 
also  depend  on  other  dimensions  and  <ru’s  dependence  on  length  need  not  be  to  the  -ft 
exponent,  or  even  particularly  close  to  this  precise  relation.  Indeed,  for  sufficiently  large 
specimens,  o  is  generally  regarded  as  a  material  properly,  independent  of  geometry  Thus 
this  simple  application  of  Griffith's  ideas  raises  a  question  as  to  whether  or  not  the 
assumption  of  a  surface  energy  term  as  the  dominant  energy  source  is  really  valid,  or  even 
sufficiently  accurate  in  certain  circumstances  The  answer,  as  far  as  cracked  geometries 
are  concerned,  has  to  lie  with  the  physical  evidence. 


The  more  modern  argument  suggested  in  support  of  the  stress  intensity  as  the 
controlling  factor  in  brittle  fracture  is  as  follows.  For  the  K— field  alone  at  the  crack-tip, 
the  associated  field  quantities  are  clearly  in  violation  of  the  underlying  yet  unpoliced 
assumptions  of  linear  elasticity,  viz,  infinitesimal  displacement  gradients,  stresses  in  the 
elastic  range,  etc.  However,  some  distance  away  from  the  crack  tip  these  K-fields  do 
comply  with  underlying  assumptions  and  may  be  regarded  as  being  valid.  To  fix  ideas,  call 
the  radius  from  the  crack  tip  outside  which  the  K-fields  are  valid  r  .  Conversely,  as  one 
approaches  the  crack  tip  for  any  configuration,  there  comes  a  point  after  which  the  K- 
fields  dominate  any  regular  fields  present  Call  this  radius  r  As  a  result,  if  r^  <  rrf,  an 

annulus  is  formed  which  may  reasonably  be  viewed  as  K  controlled,  i.e.  a  region  in  which  K 

sets  the  boundary  conditions  in  effect  for  the  process  zone  at  the  crack  tip.  Thus  under 

these  circumstances,  given  similar  fracture  processes  in  the  same  material.  K  determines 

when  fracture  occurs 

The  drawback  to  this  argument  clearly  is  that  r^  may  not  be  less  than  rrf  Some 
indication  that  this  can  occur  is  contained  in  [11]*,  [12]*,  wherein  three  very  similar  pac- 
man  geometries  under  the  same  loading  are  treated.  By  changing  the  roof  of  the  pac- 
man's  mouth  slightly.  K  switches  from  on  to  off  to  on  again.  Here  in  some  sense  rd  is 
going  from  some  finite  value,  to  zero,  to  some  finite  value.  Hence  at  some  point  in  this 
sequence,  no  K-controlled  annulus  exists.  And  further  examples  can  also  be  constructed. 
Thus  the  question  arises  as  to  whether  or  not  this  approach  is  sufficiently  reliable  to  be 
useful  in  practice:  again  the  answer  rests  with  the  physical  evidence. 

Turning  to  the  physical  evidence,  perhaps  the  most  basic  question  is  to  what  extent  is 
the  critical  value  of  K,  the  fracture  toughness,  a  material  constant  In  particular,  to  what 
degree  is  the  plane  strain  fracture  toughness,  K(c,  a  material  constant,  since  this  is  probably 
the  most  carefully  governed  by  standards  in  measurement  of  all  parameters  in  fracture 
mechanics  today  A  review  for  two  materials  which  possibly  have  had  the  most  extensive 
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testing  in  this  regard  of  all  steel  and  aluminum  alloys,  respectively,  is  given  in  [13]*  a 
summary  of  the  results  is  contained  in  Table  1  (bars  indicate  mean  values,  A  indicates  95% 


confidence  limit,  a  is  yield  stress,  K£  plane  stress  fracture  toughness).  In  sum,  K|c  is 


about  five  times  more  variable  than  the  yield  stress  and  accordingly  quite  dubious  as  a 
material  property. 


Table  1.  Critical  value  of  K  as  a  material  property 


Material 


A  <r 


Property  ranges 
AK. 


AK 


7075  T6  Al  +  11%  +81%  K,  +55%  K 

—  Y  —  Ic  —  c 

4340  Steel  +9%  +50%  K 

—  Y  —  IC 


A  second  basic  question  in  terms  of  physical  evidence  is  how  well  can  K.  be  used  to 
predict  fracture  in  the  simplest  of  situations.  Focusing  on  geometrically  similar 
configurations  subject  to  a  change  of  scale  alone  -  arguably  the  most  fundamental  of  tests 
of  predictive  capability  -  an  extensive,  if  not  comprehensive,  survey  of  all  the  pertinent 
data  is  furnished  in  1 14)*^  [15]**  .  The  answer  found  is  summarized  in  Fig  1  This 
figure  shows  the  percentage  of  times  LEFM  is  within  10%  of  predicting  fracture  -  useful 
agreement  -  and  the  percentage  of  times  it  is  within  5%  -  good  agreement  These 
percentages  are  shown  as  a  function  of  scaling  with  the  number  of  independent  tests 
involved  for  each  class  of  scale  factors  given  in  parentheses  When  scaling  is  small  (2-3) 
there  is  really  little  to  predict  even  so  K  fails  to  furnish  useful  predictions  over  70%  of 
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the  time.  When  scaling  is  moderate  (3-4),  K  fails  to  furnish  useful  results  over  80%  of 
the  time,  while  when  it  is  significant  (>4)  it  fares  somewhat  worse  until  eventually  there  are 
no  instances  of  satisfactory  performance.  And  good  predictive  capability  is  even  more 
rare.  This  performance  in  such  carefully  controlled  laboratory  situations  could  perhaps  have 
been  expected  given  the  outcome  of  the  earlier  question  concerning  K|c  as  a  material 
constant,  but  is  nonetheless  most  unsatisfactory. 
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CONCLUDING  REMARKS 

Fracture  mechanics  today  typically  is  successful  in  identifying  and  quantifying 
singularities,  prerequisites  to  a  useful  technology.  The  present  interpretation  of  the 
resulting  quantified  singularity,  however,  is  not  good  enough.  More  precisely,  the  choice  of 
the  stress  intensity  factor  as  the  parameter  governing  brittle  fracture  under  monotonic 
loading,  while  appealing  in  its  simplicity,  is  not  up  to  handling  the  complex  phenomena 
involved.  There  is  a  serious  and  immediate  need  to  improve  our  technology  here 

Further  investigations  are  warranted  regarding  the  performance  of  fracture  mechanics 
in  the  fatigue  or  cyclic  life  arena  Given  the  unsatisfactory  performance  of  fracture 
mechanics  in  the  simpler  monotonic  loading  instance,  there  would  seem  to  be  little  reason 
for  optimism  here.  Probably  significantly  different  alternatives  will  need  to  be  developed 
for  both  monotonic  and  cyclic  loading  before  engineering  can  be  provided  with  an  adequate 
technology. 
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APPENDIX  1:  ON  THE  SINGULARITIES  IN  REISSNER’S  THEORY  FOR 
THE  BENDING  OF  ELASTIC  PLATES 
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Brief  Notes 


A  Brief  Note  is  >  short  paper  that  presents  a  specific  solution  of  technical  interest  in  mechanics  but 
which  does  not  necessarily  contain  new  general  methods  or  results.  A  Brief  Note  should  not  exceed 
1500  words  or  equivalent  (\  typical  one-column  figure  or  table  is  equivalent  to  250  words;  a  one  line 
equation  to  30  words).  Brief  Notes  will  be  subject  to  the  usual  review  procedures  prior  to 
publication.  After  approval  such  Notes  will  be  published  as  soon  as  possible.  The  Notes  should  be 
submitted  to  the  Technical  Editor  of  the  Journal  or  Applied  Mechanics.  Discussions  on  the  Brief 
Notes  should  be  addressed  to  the  Editorial  Department,  ASME,  United  Engineering  Center,  345 
East  47th  Street,  New  York,  N.  Y.  10017,  or  to  the  Technical  Editor  of  the  Journal  or  ArruEO 
Mechanics.  Discussions  on  Brief  Notes  appearing  in  this  issue  will  be  accepted  until  two  months 
after  publication.  Readers  who  need  more  time  to  prepare  a  Discussion  should  request  an  extension 
of  the  deadline  from  the  Editorial  Department. 


On  the  Singularities  in  Reissner’s  Theory 
for  the  Bending  of  Elastic  Plates 

W.  s.  Barton 1  aad  G.  B.  Sinclair' 


Wedge-shaped  elastic  plates  under  bending,  with  the  edges 
forming  the  wedge  vertex  being  either  stress-free,  clamped  or 
simply  supported,  are  characterized  as  to  possible  singular 
behavior  within  the  context  of  Reissner’s  plate  theory. 

Introd  action 

Probably  the  first  singularity  analysis  of  an  angular  elastic 
plate  under  bending  is  William's  treatment  using  the  classical 
theory  ()].  In  the  classical  theory  it  is  possible  to  satisfy  stress- 
free  conditions  at  an  edge  solely  in  an  approximate  way,  since 
only  two  boundary  conditions  can  be  enforced  and  there  are 
three  stress  resultants.  As  the  boundary  conditions  play  an  im¬ 
portant  role  in  governing  singular  behavior  at  the  vertex  of 
any  comer  in  a  plate  it  is  to  be  expected  that  Reissner's  theory 
(2|,  which  admits  three,  physically-natural,  boundary  condi¬ 
tions  on  an  edge,  may  offer  an  unproved,  albeit  singular, 
representation  in  these  instances.  This  is  the  expectation  that 
possibly  motivated  other  analysts  (e.g.,  Wang  |3|)  to  perform 
analyses  of  complete,  individual,  crack  problems  using 
Reissner's  theory  rather  than  the  classical,  and  indeed  more 
physically  sensible  results  are  derived  in  these  analyses. 
Specifically,  for  the  crack-tip  on  the  tensile  face  of  the  plate 
the  same  hydrostatic  singular  field  ahead  of  the  crack  as  oc¬ 
curs  in  the  extensional  case  of  a  cracked  plate  under  tension  is 
found;  this  is  in  contrast  to  the  classical  bending  theory  in 
which  the  principal  stresses  ahead  of  the  crack  differ  in  sign 
and  magnitude.  As  a  result  it  would  seem  reasonable  to  at¬ 
tempt  the  analogue  of  Williams'  study  for  the  classical  theory 
U|  and  explore  the  singularities  in  Reissner’s  plate  theory  for  a 
wider  range  of  geometries  than  that  investigated  elsewhere  and 
for  a  full  range  of  boundary  conditions;  this  is  the  intent  of  the 
present  note. 

We  begin  by  formulating  a  class  of  problems  for  a  wedge- 
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shaped  elastic  plate  under  bending  within  the  context  of 
Reissner's  theory.  Next  we  establish  suitable  solution  forms 
for  the  dominant  asymptotic  response  near  the  wedge  vertex 
and  set  down  conditions  for  the  existence  of  these  fields.  "Die 
conditions  basically  involve  the  analysis  of  an  eigenequation 
for  each  pair  of  edge  conditions  considered.  The  note  con¬ 
cludes  by  displaying  these  eigenequations  and  discussing  the 
eigenvalues  satisfying  them  which  give  rise  to  singularities. 

Formulation 

The  plate  has  thickness  h  and  occupies  the  open  wedge 
region,  R, 

*.((r,*)K><r<«».-a/2<#<o/2|(0<a*2v),  (I) 

where  (r,t)  are  the  polar  coordinates  of  a  point  P  in  the  wedge 
with  respect  to  the  origin,  0,  at  the  wedge  vertex,  and  a  is  the 
vertex  angle  (Fig.  I). 

The  plate  is  comprised  of  a  homogeneous,  isotropic,  and 
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Table  1  Eige  nequations  for  Reissner’s  plate  theory  generating 


Edge  Conditions 

Eigenequation 

Constants 

Stress-free/stress-free 

tinXa  *  CjX 

C,  «(-)*sina 

damped  /clamped 

tinXa  «  CjX 

C}-(-  I/«)‘*'sino 

Clamped /stress-f  ree 

sin2Xa  -  Cj  -  C4X2 

Cy  *  4/e(  1  +  e)2 

C4  «(l/*)sinJo 

Stress-free /simply  supported 

sm2Xa-C5X 

Cs  •  sin2a 

Ctamped/simply  supported 

sin2Xn  =  C(X 

Ct-(-l/i)sin2a 

Simply  supponed/simply 
supported 

cosXa  *  C7 

Ci  “(- )*COSa 

Note:  km  |,j  for  symmetric  or  anti-symmetric  loading,  respectively,  and  «« (3-  »)/(!  +  r). 

linear  elastic  material  having  Young’s  modulus,  E ,  and 
Poisson's  ratio,  ».  For  this  plate  the  stress  resultants  and  rota¬ 
tions  of  Reissner’s  plate  theory,  in  the  absence  of  surface 
loading,  can  be  expressed  in  terms  of  the  out  of  plane  deflec¬ 
tion,  w(r,9),  and  a  single  stress  potential,  x(r,S).  That  is, 

y  .  1  3x  ,/  -  a* 

k“tit- K  ir¬ 


y‘; 
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_/  32w  1  9w  1  92w  \ 
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P'  D(l-»)  r  99  dr  ’ 

on  R,  where  V„  V„  the  shear  resultants,  M„  M„  Af*.  the  mo¬ 
ment  resultants  and  0(,  the  rotations,  are  functions  of  r,  9, 
defined  in  the  usual  manner,  with  y»AJ/10,  /)«£AV 
12(1  -  r2),  the  last  being  the  flexural  rigidity.  Then  the  field 
equations  of  Reissner's  theory  are  reduced  to  the  Cauchy- 
Rjemann  equations  for  the  functions  x~yV2X  and  Dv2w, 
i.e„ 

^(x-YV«*)-L-|-a>v»w), 


I  9 

T¥(*-tt,x)' 


-(DVi  w). 


on  R  Here  V2  is  the  Laplacian  operator  in  cylindrical  polar 
coordinates. 

On  each  wedge  face  three  homogeneous  boundary  condi¬ 
tions  are  to  be  satisfied.  These  three  boundary  conditions  are 
combined  in  sets  of  edge  conditions  to  model  various  edges  as 
follows: 

Stress-free  M,  *  *  0,  V,  m  0, 

Clamped  w*0,  (4) 

Simply  supported  M,  ■  0. 6, -  0,  w«  0, 

on  9  »  *  a/2  (Osr<  »).  These  three  cases  combine  to  give  six 
distinct  problems  for  the  wedge.  When  the  same  conditions 
apply  on  each  face,  it  is  possible  to  distinguish  between  tym- 
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metric  and  anti-symmetric  contributions.  Thus,  in  effect,  nine 
problems  are  considered.  Ensuring  bounded  displacements 
concludes  our  formulation  and  limits  the  singular  behavior 
admitted-while  this  formulation  is  then  still  not  complete,  it 
suffices  given  our  objective  of  characterizing  possible  singular 
fields  at  the  wedge  vertex.  We  next  consider  the  construction 
of  suitable  sets  of  asymptotic  solutions. 

Singularity1  Analysis 

From  the  form  of  the  relations  in  (2),  observe  that  if 
wwOfr**1),  x“0(ri'*l)»  as  r— 0  on  R.  X  a  constant,  then  the 
shear  resultants  and  rotations  are  0(r*),  while  the  moment 
resultants  are  Ofr*- *).  Accordingly  we  seek  to  construct 
separable  solutions  for  w,  x  satisfying  the  governing  equations 
(3)  which  furnish  six  independent  constants  multiplying  these 
dominant  contributions  as  r— 0,  thereby  providing  a  means  of 
satisfying  the  six  boundary  conditions  contained  in  any  pair  of 
the  edge  conditions  (4).  To  this  end,  and  noting  that  x  -  kv:x 
and  DV2w  are  harmonic  functions  and  the  interrelations  in 
(3),  we  therefore  take  as  our  asymptotic  solution  forms  for  x 
and  w,  the  biharmonic  functions2 

X-f^’F<k.»)  +  0<r**»), 
w«r‘‘lC(M)  +  0(r‘*,)1  (5) 

as  r—  0  on  R,  where 
FT  X.9)  “  (A,cos(X  +1)9  +  h2sin(X  +  1  )9 

+ 9jCOs(X  -  1)9  +  b,  sin(X  -  1  )9), 
C(X.9)«(h5tin(X  +  1)9  + Atcos(X+  1)9 

-yhjSin(X-  l)9  +  yb4cos(X-  1)9)//). 

The  stress  and  moment  resultants  and  the  rotations  in  (2)  may 
then  be  written  as 

F>r*F”  +  0(r*‘2),  K,--(X+l)r*F  +  0(rk*2), 
Af,«r*-'|-2yXF’-D((A+l)(l  +  rX)G  +  G*))  +  0(r*<,)1 
Af, «  r* ■ '  I2yXF'  -  Z>(  ( X  +  I XA  +  r)G  +  eC *  ))  +  (Hr* * ' ). 

Af  +  «  r*  - '  h  ( F' -  (X  +  1 KX  -  I  )F)  - />(  1  -  r)XG  ’  ]  +  (Hr*  * 1 ). 

^-^[•^—(X+DF-Gj+tK^*1),  (6) 

as  r— 0  on  R,  where  the  primes  denote  differentiation  with 
respect  to  9. 

With  this  set  of  separable  functions  for  w  and  x  the 
singularity  analysis  proceeds  in  a  manner  similar  to  that 
developed  by  Williams  (4]  for  power  singularities  and  by 
Dempsey  and  Sinclair  (3]  for  logarithmic  singularities.  Impos¬ 
ing  the  displacement  regularity  requirements  on  (6)  and  con- 

that  onr  cannot  aw  Rettsner's  >otufion  [})  and  have  a  «uffioeni 
member  of  inOcpcndtm  coasanti  available  for  the  asymptotic  analysis 
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fining  attention  to  singular  solutions,  the  resulting  conditions 
may  be  summarized  by 

M  *  (Hr** ')  for  real  X  satisfying  D  «  0, 0  <  X  <  1 , 

/  ,  f  sin(t|tor)"] , 

M-Olr*"1  ^  V)  for  complex  X»{  +  h\ 

V  Lcos(,Mr)J ' 

satisfying  0*  0, 0<  A/?eX<  1,  (7) 


M  *  (%/*  -  '/nr)  for  real  X  satisfying 

t r~mD 

- *0.  m<n,  0<Xsl, 

d\-m 


M-0  (r*-'lnr 


r  sin(i)lnr)  ~J , 

V  I  for  complex  X  *  £  +  m 
„  cos(i/nr)  J  ' 


tr-mD 


0,  mot,  Oc/IeXsl, 


for  the  first  three  cases  in  Table  I  by  Williams  [4]  and  are 
decomposed  in  effect  into  symmetric  and  anti-symmetric  parts 
by  Kalandiia  [61,  thus  accounting  for  the  first  five  cases.  The 
roots  for  the  last  case,  simply  supported/simply  supported, 
can  be  determined  by  inspection.  Some  of  the  eigenvalues 
given  in  [4,  6)  are  actually  the  real  parts  of  complex  solutions; 
however,  no  truly  comprehensive  search  for  complex  roots  ap¬ 
pears  to  be  available  in  the  literature.  Such  a  parameter  study 
is  outside  the  scope  of  the  present  work.  Nonetheless,  for  any 
given  application  the  determination  of  complex  eigenvalues 
proceeds  routinely  on  separating  the  pertinent  eigenequation 
for  the  specific  a-value  into  real  and  imaginary  parts  and  solv¬ 
ing  the  resulting,  simple,  simultaneous  pair  of  transcendental 
equations.  Likewise,  logarithmic  singularities  have  not  been 
exhaustively  searched  for,  but  are  straight  forward  to  check 
for  in  any  specific  instance. 

In  conclusion  we  remark  that  the  correspondence  between 
the  singular  fields  in  Reissner’s  theory  and  those  in  extensional 
plate  theory  is  not  restricted  merely  to  the  singular  eigen¬ 
values,  but  carries  over  to  the  actual  eigenfunctions  which 
share  the  same  r  and  8  dependences,  as  can  be  deduced  from 
the  solution  (6)  and  its  counterpart  in  Williams  [4]. 


as  r— 0  on  /?,  where  M,,  M*)  is  the  vector  of  mo¬ 

ment  resultants.  In  (7),  D  is  the  determinant  of  the  coefficient 
matrix  stemming  from  the  substitution  of  (6)  into  a  set  of  edge 
conditions  drawn  from  (4)  and  n  is  the  order  of  this  matrix,  m 
its  rank.  For  any  particular  combination  of  the  edge  condi¬ 
tions  (4)  for  a  wedge  angle  a,  the  values  of  X  in  the  ranges 
given  in  (7)  may  be  regarded  as  the  singular  eigenvalues  of  the 
eigenequation ,  0*0.  We  now  investigate  the  eigenequations 
resulting  from  such  expansions. 

Eigenequations 

For  the  particular  problem  of  the  symmetric  bending  of  a 
stress-free/stress-free  wedge,  substituting  (6)  with 
b,  •  b}  *  b,  *  0  therein  into  the  first  of  (4)  and  expanding  the 
determinant  of  the  resulting  3x3  coefficient  matrix  leads  to 

sin(X  +  I)a/2(Xsina  +  sinXa)*0(0<as2*).  (8) 

Equation  (8)  factors  into  two  equations;  however,  the  first  of 
these,  while  not  generating  a  completely  trivial  solution,  does 
not  give  rise  to  any  moment  resultants  and  therefore  con¬ 
tributes  no  singular  fields.  Consequently,  it  may  be  discarded 
leaving  as  our  eigenequation  for  this  case  only  the  second  fac¬ 
tor.  The  eigenequations  for  the  remaining  combinations  of  the 
edge  conditions  each  possess  similar,  non-singular, 
multiplicative  factors.  In  Table  I,  we  suppress  these  and  list 
only  those  pans  of  each  eigenequation  that  have  attendant 
singular  fields. 

Comparison  of  the  first  three  cases  in  Table  1  with  the  cor¬ 
responding  extensional  cases  given  in  Williams  (4)  shows  the 
eigenequations  to  be  identical.  Examining  the  founh  and  fifth 
cases  in  Table  1  and  noting  that  the  conditions  for  the  simply 
supponed  edge  in  (4)  are  the  same  as  anti-symmetry  re¬ 
quirements,  we  see  these  eigenequations  are  equivalent  to  the 
anti-symmetric  parts  of  the  first  and  second  eigenequations, 
respectively,  for  a  wedge  of  angle  2a.  It  follows  that  these  two 
cases  are  also  effectively  contained  in  Williams’  extensional 
analysis  |4).  Finally,  taking  as  the  physical  analogue  of  the 
simply  supported/simply  supported  edge  condition,  the  exten¬ 
sional  anti-symmetry  conditions,  u,*0,  o,  *0  where  u,  is  the 
radial  displacement  and  o,  is  the  tangential  normal  stress,  we 
find  that  the  last  case  too  has  a  corresponding  extensional 
eigenequation.  The  significance  of  this  correspondence  is  that 
discussions  in  the  literature  on  the  extensional  eigenequations 
are  directly  applicable  to  elastic  wedges  generated  by 
Reissner's  bending  theory. 

Solutions  for  the  dominant  singular  real  part  of  X  are  given 
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Dyuamk  Behavior  of  Beam  Structures 
Carrying  Moving  Manes 

S.  Salgal' 

latradactksa 

The  dynamic  response  of  structure*  carrying  moving  masses 
is  a  problem  of  widespread  practical  significance.  A  detailed 
survey  of  research  efforts  in  this  field  was  given  by  Stanisic  et 
al.  [2[.  The  original  problem  is  nonlinear  in  both  local  and 
convective  derivatives  [3]  and  is  complicated  by  the  presence 
of  a  Dirac-Delta  function  as  a  coefficient  in  the  differential 
equation  of  motion.  Previous  methods  [2[  applied  for  the 
solution  of  this  problem  are  approximate  in  nature  and 
tedious  in  their  hierarchy  of  mathematical  operation.  Recent¬ 
ly,  Stanisic  |3[  expressed  the  solution  in  terms  of  eigenfunc¬ 
tions  satisfying  the  boundary,  initial  and  transient  conditions, 
for  a  heavy  mass  moving  over  a  simply  supported  beam. 
However,  in  engineering  practice  there  are  problems  that  in¬ 
volve  more  complex  boundary  conditions  and,  therefore,  it  is 
of  phenomenological  interest  to  look  into  the  physics  of  the 
dynamical  behavior  of  a  clamped  and  a  cantilever  beam  under 
the  action  of  heavy  moving  masses.  The  present  study  extends 
Stanisic's  theory  [3J  to  study  the  dynamic  behavior  of  a 
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Lafayette.  Indiana  47907.  Suidcm  Member.  ASME. 
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PATH  INDEPENDENT  INTEGRALS  FOR  COMPUTING 
STRESS  INTENSITY  FACTORS  AT  SHARP  NOTCHES 
IN  ELASTIC  PLATES 

C  B  SINCLAIR.  M  OKAJIMA  AND  J  H  GRIFFIN 

Department  of  Mechanical  Engineering  and  the  Ctnier  for  the  Joining  of  Materials.  Camcgie-Mellon  University. 

Pittsburgh.  Pennsylvania.  U.S.A. 

SUMMARY 

A  set  of  path  independent  integrals  is  constructed  for  the  calculation  of  the  generalized  stress  intensity 
factors  occurring  in  elastic  plates  having  sharp  re-entrant  comers  or  notches  with  stress-free  faces  and 
subjected  to  Mode  I,  II  or  III  type  loading.  The  Mode  I  integral  is  then  demonstrated  to  enjoy  a  reasonable 
degree  of  numerical  path  independence  in  a  finite  element  analysis  of  a  test  problem  having  an  exact 
solution.  Finally,  this  integral  is  used  on  the  same  problem  in  conjunction  with  a  regularizing,  finite 
element,  procedure  or  superposition  method.  The  results  indicate  that  sufficiently  accurate  estimates  of 
these  stress  intensity  factors  for  engineering  purposes  can  be  achieved  with  little  computational  effort. 

INTRODUCTION 

In  a  number  of  engineering  structures,  sharp  re-entrant  comers  or  notches  are  introduced, 
usually  to  facilitate  fabrication.  At  the  vertices  of  these  notches,  stress  concentrations  arise 
making  them  likely  sites  for  fatigue  crack  initiation  and  therefore  the  potential  sources  of 
ultimate  failure.  Elastic  analyses  of  such  configurations  result  in  stress  singularities  (see,  e.g.. 
References  1  and  2).  While  these  singularities  are  physically  unrealistic  in  themselves,  it  is 
possible  that,  with  a  sufficiently  accurate  determination  of  their  nature  and  participation,  an 
extension  of  the  now  accepted  notion  of  a  stress  intensity  factor  X  might  prove  of  value. 
Although  such  a  generalization  lacks  the  thermodynamic  argument  in  terms  of  energy  release 
rates  that  underlies  the  physical  significance  of  X  for  cracks,  it  could  lead  to  an  improvement 
over  present  practice  in  initiation  calculations  for  notches.  One  current  approach  to  such 
calculations  is  to  assume  that  a  crack  has  actually  formed  at  the  notch  tip,  then  compute  the 
usual  stress  intensity  factor  variation  under  cyclic  loading,  AX,  and  compare  it  with  material 
threshold  values  to  see  if  crack  growth  commences:  with  a  generalized  X  for  the  notch  itse'f, 
the  AX  could  be  computed  directly  and  compared  with  an  accompanying  set  of  experimental 
threshold  values  for  the  notch.  In  this  way,  increased  life  could  be  detected  in  circumstances 
in  which  load  levels  are  too  low  for  the  formation  of  a  crack  at  the  notch  tip.  The  intent  of 
this  investigation  is  to  contribute  to  the  analytical  component  of  such  a  technology  by  developing 
an  efficient  computational  procedure  for  accurately  assessing  generalized  stress  intensity  factors 
at  stress-free  notches  in  elastic  plates,  one  which  is  readily  applicable  to  the  wide  range  of 
complex  configurations  encountered  in  practice. 

The  basis  of  our  approach  to  meet  the  foregoing  objective  lies  in  the  development  of  a  set 
of  path  independent  integrals  which  pick  of!  the  X  at  notches.  Such  path  independent  integrals 
for  singular  elasticity  problems  probably  owe  their  origin  to  Eshelby's  work1  in  the  early  1950s, 
and  a  number  of  different  forms  for  these  integrals  have  evolved  over  the  intervening  years 
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(e  g.  Sanders'  integral;4  see  also  Freund*  and  the  references  contained  therein).  All  of  the 
foregoing  integrals  apply  to  a  crack.  Recently  Stern  and  Sonifc  developed  an  integral  for  a 
right -angle  corner  under  fixed-free  conditions;  here  we  follow  the  ideas  in  Reference  6  and 
construct  a  set  of  like  integrals,  termed  H  integrals,  for  a  stress-free  notch  in  an  elastic  plate. 

The  computational  advantage  of  path  independent  integrals  in  the  analysis  of  singular 
elasticity  problems  stems  from  the  fact  that  the  only  numerical  errors  accrued  in  their  calculation 
derive  from  the  numerical  approximations  being  used — there  is  no  additional  source  of  errors 
such  as  that  due  to  the  truncation  of  an  eigenfunction  expansion  as  in  local  fitting  methods.  It 
is  to  be  expected  therefore  that,  when  used  in  conjunction  with  a  finite  element  method  which 
recognizes  and  handles  the  singularity  present  or  some  equivalent  numerical  approach,  a 
procedure  results  which  is  easily  adapted  to  varying  configurations  yet  which  yields  the  necessary 
resolution  for  practical  purposes  in  return  for  a  modest  computational  effort.  Such  is  indeed 
demonstrated  to  be  the  case  in  the  second  section  of  the  paper,  wherein  the  numerical  path 
independence  of  an  H  integral  drawn  from  the  previous  section  is  examined,  following  which 
the  same  integral  is  used  in  conjunction  with  the  regularizing  procedure  of  Sinclair  and  Mullan.7 
The  paper  then  closes  with  some  remarks  on  other  applications  and  extensions. 

CONSTRUCTION  OF  THE  PATH  INDEPENDENT  INTEGRALS 

Here  we  first  focus  on  the  in-plane  symmetric  loading  of  a  notch  in  an  elastic  plate  and  develop 
a  path  independent  integral  for  the  stress  intensity  factors  in  this  class  of  problems.  We  then 
discuss  other  integrals  for  anti-symmetric  and  out-of-plane  configurations. 

To  formulate  our  initial  class  of  symmetric,  in-plane,  notch  problems  we  consider  a  plane 
finite  wedge  subtending  an  angle  2 a  at  its  vertex  ( n/2<a<ir ),  let  (x,,x2)  be  rectangular 
cartesian  co-ordinates  having  origin  at  the  wedge  vertex,  the  x,-axis  bisecting  the  angle  there, 
and  take  R  to  be  the  open  region  defined  by  that  part  of  the  wedge  in  the  upper  half-plane, 
with  aR  being  its  boundary  (Figure  1).  To  facilitate  further  the  formulation  and  subsequent 
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analysis,  we  introduce  companion  cylindrical  polar  co-ordinates  (r,  6)  related  to  (x,,  x2)  by  the 
transformations, 

x,  =  r  cos  6  x2  =  rsin  0  (Os  r<oo,  -  rr  <  0s  it)  ( 1 ) 

Then,  in  general,  we  seek  the  stresses  <r„  and  displacements  «,(/,;=  1,2)  throughout  R  as 
functions  of  x,.  x2  satisfying  the  following.  The  two-dimensional  stress  equations  of  equilibrium 
in  the  absence  of  body  forces, t 

<t„i  =  0  on/?  (2) 

The  stress-displacement  relations  for  a  homogeneous  and  isotropic,  linear  elastic  plate, 

=  m[((3- n)/(k-  Dlfi.jUnn  +  uM+  u,,]  on  R  (3) 

wherein  n  is  the  shear  modulus.  *  equals  3-4v  for  plane  strain  and  (3-v)/(l  +  v)  for  plane 
stress  t-  being  Poisson's  ratio,  and  6,,  is  the  Kronecker  delta.  The  stress-free  conditions  on  the 
notch  face. 

0  on  HR  when  x,/x2  =  cot  a  (4) 

where  n,  are  the  components  of  the  unit  outward  normal  to  A/?.  The  symmetry  requirements 
ahead  of  the  notch  tip. 

ct,2  =  0  u2  =  0  on  A/?  when  x2  =  0  (5) 

The  boundary  conditions  prescribing  the  tractions  or  displacements  on  the  remainder  of  A/?, 

a1(n(  =  s?  on  A,/?  u,  =  u?  on  f>2R  (6) 

wherein  sf,  u?  are  the  applied  tractions  and  displacements  and  A,/?,  a2/?  are  complementary 
subsets  of  A/?  excluding  the  intervals  in  (4),  (5).  And  finally  the  regularity  requirements  at  the 
notch-tip  which  insist  that  the  displacements  are  bounded  there, 

«,*0(1)  on/?asr-»0  (7) 

Specifically  we  seek  to  extract  from  the  stress  and  displacement  fields  meeting  these  require¬ 
ments  the  generalized  stress  intensity  factor  present  defined  by 

K,  m  lim.  V(2ir)xJ"*o22  onx2*0  (8) 

*,-•0 

Here  the  subscript  I  denotes  the  opening  symmetric  mode.  Mode  1,  and  A  is  the  singular 
eigenvalue  characterizing  the  only  singular  stress  field  possible  at  the  re-entrant  comer  as 
identified  after  Williams.  :t  That  is,  A  is  the  root  of  the  transcendental  eigenequation, 

sin  2Aa  *  —A  sin  2a  (0  <  A  <  1 )  (9) 

with  an  associated  eigenfunction  having  stresses  and  displacements  in  the  neighbourhood  of 
the  notch  tip  behaving  in  accordance  with 

o,, « Ofr*'1)  «, -0(rA)  on/?8sr-*0  (10) 


r  Although  the  lowet  cue  lubicnptt  only  range  over  the  integers  (1,2),  the  usual  index  notation  conventions  apply: 
repeated  subscripts  imply  summation  and  a  subscript  preceded  by  a  comma  indicates  partial  differentiation  with 
respect  to  the  corresponding  cartesian  co-ordinate 
t  See  Gregory*  lor  a  completeness  argument  lor  Williams'  eigenlunctions. 
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To  construct  a  path  independent  integral  for  computing  Kx  we  proceed  as  follows.  In  general, 
path  independent  integrals  can  be  devised  by  invoking  the  divergence  theorem  of  Gauss  on 
ensuring  the  divergence  of  the  integrand  is  zero:  in  elasticity,  this  can  be  done  in  effect  with 
Betti’s  reciprocal  theorem  (refer,  e.g.,  to  Reference  9,  p.  355)  which,  in  a  plane,  can  be  stated 

|  (<r„u*  -<r*u,)n,ds  =  0  (11) 

In  (11),  the  integration  is  performed  in  a  counter-clockwise  sense  around  T  which  is  any  closed 
contour  in  R,  a J,  u?  are  complementary  fields  satisfying  the  same  field  equations  as  a,,,  u„ 
namely  (2),  (3),  n,  are  now  the  components  of  the  unit  outward  normal  to  T,  and  ds  is  an 
infinitesimal  line  segment  of  T.  Next  we  choose  T  as  the  contour  comprised  of  any  path  X 
which  emanates  on  the  lower  notch  flank  and  terminates  on  the  upper,  a  circular  path  of  radius 
r  from  the  upper  flank  to  the  lower,  and  two  closing  straight  paths  along  the  flanks  (Figure 
2)  On  letting  the  radius  of  the  inner  circular  arc  shrink  to  zero  it  is  clear  that  only  those  parts 
of  the  integrand  which  behave  like  0(l/r)  as  r-*0  can  contribute  to  this  portion  of  the  integral 
in  (11).  In  order  that  such  contributions  stem  from  the  singular  K,  fields  alone  we  therefore 
require  that 

er*  ■  0(r~*_l),  u*-0(r'*)  oni?asr-*0  (12) 

Then  (10),  (12)  imply  that  the  product  of  the  Kt  terms  with  the  complementary  field  in  (11) 
has  the  desired  1/r  behaviour  while  all  other  fields  in  a,r  u,  contribute  terms  o(l/r)  as  r-» 0 
by  virtue  of  the  fact  that  the  singularity  is  the  dominant  one  present.  Unfortunately,  under 
this  limiting  procedure  the  possibility  of  divergent  integrals  along  the  flanks  arises  because  of 
the  potential  of  1/r  terms  there.  Fortunately,  if  A  is  an  eigenvalue  satisfying  (9)  so  is  -  A,  so 
that  the  complementary  fields  can  be  further  required  to  be  eigenfunctions  themselves,  thereby 
rendering  both  a,,n,  and  <r*n,  zero  on  the  flanks  and  making  the  integrals  there  identically 
zero.  Under  these  conditions  the  integral  around  X  equals  the  counter-clockwise  integral  around 
the  circular  arc  in  the  limit  as  r-*0  which,  with  the  selection  of  appropriate  participation  factor 
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in  the  complementary  eigenfunction,  can  be  adjusted  to  recover  X,.  That  is,  we  define  //,  by 


Ht= ~a*u‘)n> ds 


and  scale  the  complementary  eigenfunction  satisfying  (2)-(5),  (12)  so  that 

(14) 

The  actual  forms  for  the  complementary  eigenfunction  are  most  readily  established  in  polar 
co-ordinates  and  can  be  derived  along  the  lines  of  Williams'  analysis2  yielding 

tr*  =-Xrr'‘_,[(A+  3)  cos  (A  +  1)0-0  cos  (A  -1)6] 

<rj*  =  K*r~A_l[(A  - 1)  cos  (A  + 1)0  -  /3  cos  (A  - 1)9] 

«-Kfr~*''[(A  +  l)sin(A  +  l)0-0  sin  (A -1)0]  (15) 

K*r~A 

u ?  *=  [(A  +  k)  COS  (a  +  1)0-0  COS  (A  -  1)0] 

2mA 

u*  «  — — —  [(A  -*)  sin  (A  +  1)0-0  sin  (A  —  1)0] 


V(8jr)>iA(A  sin2  a  +  sin2  Aa) 


(1  +  w)(sm  2a  +  2a  cos  2 Aa) 


A  cos  2a  +cos  2Aa 


Before  moving  on  to  the  application  of  the  H,  integral  of  (13),  several  remarks  concerning 
its  use  and  some  simple  extensions  are  in  order.  First,  because  of  the  symmetry  involved,  it  is 
clear  that  H,  can  be  computed  as  twice  the  integral  around  that  part  of  I  within  R  alone. 
Second,  (14)  holds  in  the  limit  a-*  it  and  the  notch  becomes  a  crack.  Third,  the  pure  traction 
conditions  on  a,/?  and  the  pure  displacement  conditions  on  B2R  are  solely  for  convenience  in 
the  formulation  and  can  be  relaxed  to  accommodate  admissible  mixtures  of  traction  and 
displacement  components.  Fourth,  extension  to  the  antisymmetric  or  Mode  11  case  is  straightfor¬ 
ward  and  leads  to 


(<r„u*- **“.)«, 


where  Ku  is  the  Mode  11,  stress  intensity  factor  defined  as  the  natural  generalization  to  a  notch 
of  the  Mode  II  factor  for  a  crack,  the  a,,,  u,  components  now  fulfil  antisymmetry  requirements 
instead  of  symmetry,  i.e. 

<7>2“0  «i“0  on  BR  when  x2  m 0  (17) 

and,  in  polar  co-ordinates,  the  antisymmetric  complementary  fields  are  given  by: 

<r*  sin  (A  — 1)0  — (A  +  3)  sin  (A  + 1)0] 

a*„  -  Kiir~k~'[p'  sin  (A  — 1)0- (A  —  1)  sin  (A  +  1)0] 

<rJi  m  K*ir~A~'[P'  cos  (A  —  1)0  - (A +  1)  cos  (A  + 1)0]  (18) 

K*r~A 

u*  »  "  -  [/?'  sin  (A  - 1)0- (A  +  «)  sin  (A  +  1)0] 


u.*“  -- 


[P'  COS  (A  —  1)0 -(A  —  k)  cos  (A  +  1)0] 


rww7:r/.7.vv.v.v..>.v.  *■  -.v 


I 


b 
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with 


Kf, 


y/(8n)tiA(A  sin5  o  -sin5  Ao) 


(1  +  *)(sin  2a -2a  cos  2Aa) 


0' 


a5-i 


A  cos  2o  -  cos  2  At* 


where  A  is  now  the  root  of  the  antisymmetric  eigenvalue  equation,  namely  (9)  with  the  minus 
sign  changed  to  plus.  Fifth,  for  problems  entailing  a  mixture  of  Modes  I  and  11,  the  fields  need 
to  be  separated  into  their  respective  symmetric  and  antisymmetric  parts  in  accordance  with 


i<r=[«1(xJ.xi)=F(-)iu1(x1.-x2)]/2  (19) 

where  u *  and  the  stresses  derived  from  substituting  u *  into  (3)  constitute  the  symmetric  part, 
etc.,  the  antisymmetric,  then  treated  individually  using  (13)  and  (16),  respectively.  Sixth 
and  last,  extension  to  include  the  outstanding  out-of-plane  mode.  Mode  III,  also  proceeds 
routinely  and  results  in 


(<T}jU)  ds; 


(20) 


where  Ktu  is  the  generalized.  Mode  Ill,  stress  intensity  factor  for  a  notch,  the  a3/,  u,  components 
are  now  the  out-of-plane  responses  (/»  1,2),  and,  in  polar  co-ordinates  (x3- r),  the  out-of¬ 
plane  complementary  fields  are: 


‘  in 


Hit  , 

2 arU‘  "V(2sr)r'*'/5* 


sin 


ird 


2a 


Kfu  pit2 

■ - cos —  Kin  * - 

■J(2ir)rl*m/2“  2a  2o5 


(21) 


TEST  PROBLEM 


In  this  section  a  plane  elastic  notch  problem  having  an  exact  solution  is  set  up  and  analysed 
using  a  finite  element  approach  and  the  integral.  The  results  found  enable  the  numerical 
path  independence  of  the  integral  to  be  examined  and  its  efficiency  as  a  computational  device 
assessed. 

For  simplicity  we  fix  attention  on  a  problem  drawn  from  the  restricted  class  formulated  in 
the  previous  section.  One  means  of  setting  up  a  problem  within  this  class  having  an  exact 
solution  is  to  superimpose  symmetric  eigenfunctions  for  a  notch  on  some  region  R,  and  we 
adopt  this  simple  device  here.  In  order  to  thoroughly  evaluate  the  performance  of  the  integral 
when  used  with  the  simplest  of  finite  elements,  the  constant -stress  triangle,  we  select  for 
combination  the  first  three  eigenfunctions  for  the  notch  since  these  contain:  the  necessary 
singular  stress  field;  a  stress  field  which,  while  continuous,  is  not  continuously  differentiable 
and  may  therefore  be  regarded  as  just  barely  being  a  regular  field  for  constant-stress  elements; 
and,  lastly,  a  stress  field  which  is  relatively  smooth  and  is  continuously  differentiable. t  More 
precisely,  we  consider  a  plate  with  a  90*  re-entrant  corner  (n  ■  3w/4)  and  take  R  to  be  the 
trapezium  within  the  upper  half-plate  (Figure  3)  defined  by 


R  -  {(x„  x2)j  -  x2  <  x,  <  a,  0  <  x2  <  a) 


(22) 


t  By  a  ' regular '  field  with  respect  to  a  given  element  we  mean  one  which  does  not  impair  the  normal  maximum 
convergence  the  element  can  enjoy;  thus  fields  are  regular  lor  analysis  using  a  constant  -stress  element  H  the  stresses 
are  continuous  (tee  Reference  10,  Section  2.2,  for  a  full  discussion  of  the  regularity  requirements  for  various  finite 
elements). 
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Figure  3.  Geometry  for  the  ten  problem,  coaru  grid  and  integration  paths 


On  this  region  we  are  to  place  the  first  three  symmetric  eigenfunctions  for  the  90s  notch, 
admissible  in  the  sense  of  (7),  and  associated  with  the  first  three  eigenvalues  of  (9)  with 
a  «  3ir/4  there  and  having  Re  A  >0.  The  forms  for  superposition  are  therefore  those  of  (15) 
with  K  J  therein  being  replaced  by  the  constants  K,,  C2,  C3  and  A  being  exchanged  for  -  A 
then  set  equal  to  the  three  corresponding  eigenvalues.  A,,  A2,  A3,  provided  the  real  parts  of 
the  expressions  in  (15)  are  taken  in  the  event  of  a  complex  eigenvalue.  In  the  interests  ol 
brevity  we  suppress  details  of  these  fields  here  and  merely  list  the  required  eigenvalues 
determined  numerically  from  (9): 

A,  -0  544  Aj- 1-629 +  i0-231  A,- 2-972 +i  0-374  (23) 

The  corresponding  eigenfunctions  satisfy  the  two-dimensional  field  equations  of  elasticity  (2), 
(3),  the  stress-free  boundary  conditions  on  the  upper  notch  face  (4),  and  the  symmetry 
conditions  ahead  of  the  crack  (5).  In  combining  them  we  first  adjust  the  participation  of  the 
A  ,-fields  so  that 

K,/<rV(2ir)a'-*'«l  (24) 

Here  X,  is  the  generalized  stress  intensity  factor  defined  in  (8)  and  <r°  equals  <r,,  for  A,  at 
x,  0.  Then,  to  complete  our  test  problem,  we  add  together  all  three  eigenfunctions 

with  the  Aj-,  A, -fields  having  participation  factors  C2,  C3  such  that  they  share  the  same  <rM 
magnitude  at  x,«a,  x2«0  as  the  A, -field  (namely  <r°)  and  take,  as  the  boundary  conditions 
on  the  remainder  of  BR  (x,  ~a,0<  x2<  a;  -  a  <x,  <a,x2«o),  the  tractions  that  result  there 
from  their  superposition.  In  this  manner  all  three  fields  contribute  to  a  comparable  degree  to 
the  test  problem. 

The  application  of  a  constant-stress  triangle,  finite  element  code  and  the  calculation  of  the 
Ht  integral  for  this  test  problem  proceeds  routinely.  Since  the  issue  of  interest  is  the  computa¬ 
tional  performance  of  the  integral,  we  discretize  R  uniformly  with  a  set  of  45*  isosceles  triangles. 
To  examine  convergence  we  employ  a  sequence  of  three  grids — a  coarse,  a  medium  and  a 
fine — with  the  medium  and  fine  grids  being  formed  from  the  coarse  by  successively  halving 
the  element  sides.  The  coarse  grid  has  48  elements  and  35  nodes  (Figure  3),  the  medium  192 
elements  and  1 1 7  nodes,  and  the  fine  768  elements  and  425  nodes.  To  examine  path  inde¬ 
pendence  we  use  three  contours  I„(n  ■  1, 2, 3)  which  have  nodes  common  to  all  three  grids 
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(Figure  3)  and  along  these  contours  use  the  mid-point  rule  to  evaluate  the  integral  in  (13) — a 
quadrature  rule  of  like  accuracy  to  that  of  the  finite  element  discretization.  The  results  so 
found  are  presented  in  the  first  half  of  Table  I. 


Table  I.  Numerical  values  of  the  dimensionless  stress  intensity  factor  for  the  test  problem  (exact  value  1 ) 


Direct  calculation  using  (13)  Calculation  using  regularizing 

on  three  integration  paths  procedure  (25)  on  three  paths 

Grid  X]  2-2  Xi  Ij  I] 


Coarse 

1044 

0-843 

0-883 

0-979 

0-975 

0  978 

Medium 

0  913 

0-939 

0-940 

0-992 

0  994 

0  994 

Fine 

0972 

0-973 

0972 

0-998 

0998 

0  999 

Of  course  the  problem  at  hand  is  tingular  and  accordingly  results  from  the  foregoing  direct 
treatment  can  be  expected  to  be  poor;  needed  is  an  approach  which  recognizes  the  singularity 
piresent  and  takes  it  into  account  if  better  results  are  to  be  realized.  One  procedure  with  this 
attribute  is  described  by  Sinclair  and  Mullan7  and  belongs  to  a  class  of  superposition  methods 
which,  in  essence,  treat  the  singular  part  of  any  problem  analytically  with  scant  regard  for 
boundary  conditions  remote  from  the  singularity  source,  then  use  a  simple  numerical  method 
to  complete  the  satisfaction  of  the  boundary  conditions  in  what  is  now  a  regular  problem. 
What  basically  distinguishes  the  procedure  in  Reference  7  from  other  superposition  methods 
is  the  use  of  path  independent  integrals  to  balance  the  analytical  and  numerical  contributions 
and  it  thus  seems  natural  to  apply  it  here.  The  end  result  so  far  as  the  stress  intensity  factor 
is  concerned  is  (see  Reference  7  for  details  of  the  surrounding  argument) 

K,*  A//*.’  (25) 

wherein  H,  is  the  numerical  estimate  of  Ht  for  the  problem  of  interest,  /?,'  the  numerical 
estimate  of  HI  for  an  associated  singular  problem  formed  by  taking  corresponding  boundary 
conditions  on  9R  as  prescribed  in  effect  by  the  singular  eigenfunction  alone  with  a  participation 
factor  of  unity.  What  (25)  says,  in  effect,  is  not  that  either  Hi  or  Rl  represent  accurate  estimates 
of  Hi  or  HI,  respectively,  but  rather  that  they  constitute  roughly  equally  crude  estimates  with 
the  errors  cancelling  in  large  part.  Central  to  this  cancellation  is  the  use  of  an  underlying 
method  for  calculating  K,  which  has  predictable  errors  due  solely  to  the  discretization — hence 
the  need  for  the  path  independent  integrals  in  (25).t 
Application  of  (25)  to  the  test  problem  is  straightforward:  the  R,  values  that  have  already 
been  computed  are  simply  divided  by  the  numerical  values  of  Ht  found  on  the  same  grid  and 
path  for  the  singular  problem  in  which  the  boundary  conditions  on  9R  are  established  by  the 
placing  of  the  A, -fields  alone  on  R  with  a  dimensionless  participation  factor  of  unity,  i.e.  as 
in  (24).  The  results  so  found  are  presented  in  the  second  half  of  Table  I 
The  numbers  in  Table  1  illustrate  the  degree  to  which  the  Ht  integral  is  numerically  path 
independent  The  variation  between  the  different  paths  is  of  the  order  erf  20  per  cent  for  the 
direct  coarse  analysis  and  converges  to  approximately  0- 1  per  cent  for  the  fine  grid;  the  variation 
between  paths  is  yet  smaller  for  the  regularizing  procedure  (25)  with  the  differences  being 
actually  less  than  0-5  per  cent  for  the  coarse  grid,  0-04  per  cent  for  the  fine.  Moreover  these 


t  Set  Reference  7  lor  t  fuller  dweuwon  of  the  advantages  of  path  independent  integrals  relative  to  netbods  aocti  at 
local  collocation  in  this  regard. 
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results  are  the  worst  experienced  in  a  number  of  numerical  experiments  since  the  signs  of  the 
participation  factors  here  are  such  that  errors  accrue  rather  than  compensate.  Accordingly  it 
would  appear  that  the  choice  of  integration  path  is  not  a  major  concern  in  applying  the  H, 
integral,  especially  when  doing  so  in  conjunction  with  the  regularizing  procedure.  The  numbers 
in  Table  1  also  indicate  the  sort  of  convergence  to  be  expected  with  H,  integral.  To  illustrate 
this  more  dearly  we  model  the  errors  for  each  analysis  with 

e  -  e0h‘  (26) 

where  e  is  the  absolute  percentage  error  in  the  dimensionless  stress  intensity  factor  for  a 
dimensionless  grid  size  of  h  :  e0  is  thus  the  error  for  unit  grid  size  and  reflects  the  initial  accuracy 
of  an  analysis,  while  the  constant  c  gauges  its  rate  of  convergence.  We  fit  (26)  to  the  medium 
and  fine  grid  results  in  Table  I  since  we  expect  that  c  is  more  likely  to  have  converged  to  a 
constant  itself  for  these  grids.  To  check  (26)  as  an  error  model  we  then  use  it  to  predict  the 
errors  for  the  coarse  grid;  with  the  exception  of  the  direct  calculation  of  H,  on  £,  these 
predictions  agree  well  with  the  actual  errors  supporting  the  use  of  (26).  The  results  of  applying 
it  can  be  summarized  by  the  average  values  determined.  For  the  direct  calculations, 


e0-13(%)  c-11 

(27) 

For  the  regularizing  procedure  results. 

e0-3<%)  c-21 

(28) 

Equations  (27),  (28)  show  that  for  this  test  problem  the  regularizing  procedure  not  only  leads 
to  markedly  better  convergence  as  expected,  but  in  addition  is  about  a  factor  of  four  more 
accurate.  Further,  the  regularizing  procedure  was  observed  to  enjoy  a  similar  superiority  over 
direct  computation  in  other  test  problems  and  the  sort  of  relative  improved  performance 
demonstrated  in  (27),  (28)  is  characteristic  of  that  obtained  in  earlier  applications  of  the 
procedure  to  crack  problems  (refer  to  Reference  7  wherein  average  values  of  e0  and  c  are  42 
per  cent  and  0-15  for  a  direct  calculation  using  the  J  integral,  in  contrast  to  8  per  cent  and 
1-4  for  the  procedure).  It  would  thus  seem  that  the  use  of  the  path  independent  integrals 
developed  here,  together  with  a  regularizing  procedure  as  in  (25),  offers  an  attractive  and 
efficient  computational  approach  for  singular  elastic  notch  problems. 

CONCLUDING  REMARKS 

The  path  independent  integrals  described  furnish  engineers  with  a  useful  tool  for  determining 
generalized  stress  intensity  factors  in  elastic  notch  problems.  This  is  especially  so  when  used 
I  with  a  regularizing  procedure  such  as  that  of  Reference  7;  then  accurate  results  can  readily 

be  obtained  using  any  standard,  constant-stress,  finite  element  code.  Improved  usage  over  that 
demonstrated  here  can  be  brought  about  by  introducing  mesh  refinement  with  the  grid  gradation 
being  tuned  to  the  complete  problem  in  the  instance  of  a  direct  calculation,  and  to  the  residual 
regular  problem  when  employing  the  procedure  in  Reference  7. 

Extension  of  these  integrals  to  treat  problems  involving  loaded  notch  flanks  is  elementary 
I  via  superposition  when  such  loads  are  constant  but  otherwise  awkward.  Extension  to  composite 

configurations  is  possible  at  the  expense  of  some  algebra  and  is  being  undertaken  at  this  time; 
extension  to  anisotropic  notch  problems  is  possible  in  principle  but  does  require  significant 
algebra.  Extension  to  three-dimensional  geometries  in  which  the  notch  tip  traces  out  a 
sufficiently  smooth  curve  is  also  possible  in  the  light  of  Aksentian’s  analysis"  which  shows  the 
singular  character  in  such  instances  to  be  that  of  the  two-dimensional  problem.  However,  for 
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geometries  that  are  more  significantly  three-dimensional,  such  as  the  intersection  of  a  notch 
with  a  free  surface,  the  integrals  given  here  would  not  be  appropriate  and  the  development 
of  suitable  new  integrals  would  appear  to  represent  a  considerable  analytical  task. 
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A  REMARK  ON  THE  DETERMINATION  OF  MODE  I  AND  MODE  II  STRESS  INTENSITY 
FACTORS  FOR  SHARP  RE-ENTRANT  CORNERS 
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The  Independent  development  in  Carpenter  (1]  and  Sinclair,  Okajima, 
and  Griffin  [2]  of  path  Independent  integrals  for  the  stress  Intensity 
factors  at  sharp  re-entrant  corners  provides  fracture  mechanics  with  a 
useful  analytical  tool.  Prompted  by  Carpenter's  recent  sequel  (3]  to  [1], 
in  which  he  extends  his  earlier  work  so  as  to  be  able  to  individually 
determine  mode  I  and  II  factors,  the  purpose  of  this  note  is  to  point  out 
that  such  a  capability  is  already  available  in  the  Integrals  of  (2).  This 
attribute  may  well  be  obscured  in  (2]  by  the  overstatement  there  that 
mixed  mode  problems  "need"  to  be  split  into  their  symmetric  and  anti¬ 
symmetric  parts  prior  to  applying  the  corresponding  path  Independent  in¬ 
tegral.  While  this  is  certainly  one  approach  for  distinguishing  between 
mode  I  and  mode  II,  it  la  by  no  means  necessary.  Simply  applying  the 
Integral  for  mode  1  in  [2]  to  a  mixed  problem  yields  the  mode  I  intensity 
and  vice  versa.  This  is  basically  because  the  complementary  fields  in 
the  mode  I  integral  are  symnetrlc  so  that  when  they  are  multiplied  by  any 
antisymmetric  fields,  including  the  singular  mode  II  field,  then  inte¬ 
grated,  no  contribution  results.  The  converse  holds  true  for  the  mode  II 
Integral.  In  detail  the  argument  is  as  follows. 

Consider  an  elastic  plate  with  a  re-entrant  corner,  or  sharp  notch, 
which  has  stress-free  flanks  (Fig.  1).  Let  x.,  x_  be  rectangular 
cartesian  coordinates  aligned  such  that  the  negative  x.-axis  bisects  the 
opening  angle  at  the  corner  into  two  angles  of  magnitude  *  -  a.  For  this 
geometry  we  can  define  the  respective  stress  intensity  factors  for  node  I, 
mode  II  by 

Rj  -  U"  '(.2')  *11"X  °22,  Kjj  •  Li**  /(2»)  x*'*  o ^  (1) 


on  X.  “  0.  In  (1),  o 
the  Stress  field  o 
values  steaming  from  their  corresponding  elgaoequatlons, 

sin  2l+o  *  -  l^sin  2®,  sin  2*  ®  ■  X  sin  2®  (2) 

with  0  <  *•*,  1”  *  1.  The  elastic  fields  associated  with  1*  and  1~  are, 
respectively,  symmetric  and  antisymmetric  about  the  x  -a^ls.  That  lsa  if 
o . . ,  u  are  the  stress  and  displacement  components  fof  1  with  o.T,  u. 
being  those  for  *  (  2 


is  the  normal  component_ln  the  x_-directlon  of 
fl,  j  “  1,  2),  while  A*,  are  the2aingular  elgen- 
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"  *2^  “  *2)*  u^(x1,-x2)  -  -  (-)1  u}(xr  x2)  (3) 


and 


°ij (xl*  ”  *2*  OjjCxj.  *2>.  u^Cxj,  ~  x2>  -  (-)1  «^(x1,  x2) 


(4) 


Now  from  f2]  we  have,  as  our  path  independent  integrals  for  determining 
the  stress  intensity  factors  for  these  singular  fields, 


HI  “  '  <0ij  UI  *  °ij  uj)njds 


(5) 


and 


Hn  •  [  <oij  \  -  °Z  uj)njds 


(6) 


Here,  £  Is  any  contour  within  the  plate  commencing  on  the  lower  notch 
flank  and  terminating  on  the  upper;  n  are  components  of  the  unit  outward 
normal  to  £  and  ds  Is  an  Infinitesimal  element  of  Its  length;  (he  Inte¬ 
gration  Is  to  be  performed  In  a  counterclockwise  sense;  and  o..,  u.  are 
the  stresses  andAjisp^|cements  of  the  symnetric  complementary  singular 
eigenfunction,  o  ,  u  those  of  the  antisymmetric  (see  Appendix  for 
details).  The  participation  factors  of  these  last  are  adjusted  in  [2] 
so  that  they  pick  off  the  corresponding  stress  intensity  factor,  l.e.. 


“l  -  V  Hn  -  *11 


(7) 


Turning  to  their  application  in  mixed  mode  problems  which  have  both 
symnetric  and  antisymmetric  singular  fields  present,  we  first  focus  on 
the  use  of  Hj  established  for  any  problem's  elastic  fields  in  [2],  the 
general  contour  £  can  be  exchanged  for  the  here  more  convenient  path  T 
(Fig.  1).  Applying  H_  of  (5)  on  T,  then  changing  variables  so  that  seg¬ 
ments  in  the  lower  half-plane  (x  <  0)  are  expressed  in  terms  of  those  in 
the  upper,  we  have 


*2  —  «  •  _  *  - 
Hj  -  r  I(o11(x1,  x2)uJ(x1,  x2)  +  °11(x1,  ”*2>ui(xl *  *  *2^ 


*  - 


*  . 


”<°il <xi *  *2>u1(xi»  *2>  +  aii<xi*  “  x2)u1(x1,-*2»)dx2 

(8) 


♦Although  the  lower  case  subscripts  only  range  over  the  Integers  (1,  2)  the 
usual  summation  convention  for  repeated  subscripts  still  applies  over  this 
reduced  range. 
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+  I  'l<012<*1.  x  )u  (x  .  52)  -  oi2(xi»  ~  *2^ui^xi‘  "  x2^ 
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-  x2^ui^xl '  x2^  ”  °i2  <xi»  ~  x  )u  (x  -  x2))]dx1. 

Now  introducing  into  (8)  the  antisymmetry  relations  satisfied  by  the  anti¬ 
symmetric  parts  of  o  ,  u  ,  namely  the  equivalent  of  (4),  in  conjunction 
with  the  symmetry  reiationg  satisfied  by  the  symmetric  complementary 
singular  eigenfunctions  o.,,  u  ,  vi2.,  the  equivalent  of  (3),  we  see  that 
all  the  combinations  in  parentheses  ()  cancel.  It  follows  that  the  only 
nontrivial  contributions  to  H  must  come  from  the  symmetric  parts  of  o  . , 
u  .  Further,  as  proven  in  [2j,  only  the  singular  part  of  these  symmetric 
fields  actually  contributes,  recovering  the  first  of  (7).  Analogously, 
may  be  shown  to  pick  off  only  Kjj  in  mixed  mode  problems. 

Thus  for  general  problems  containing  a  mixture  of  sysnetrlc  and  anti¬ 
symmetric  fields,  in  addition  to  the  strategy  of  separating  the  given 
problem  into  its  symmetric  and  antisymmetric  constituents  if  it  occupies 
a  region  geometrically  symmetrical  about  the  x^-axls  ,  one  can  simply 
apply  the  H..  and  H_..  Integrals  directly  on  the  path  of  choice  to  discern 
K_  and  respectively.  Either  technique  provides  reliable  estimates  of 
the  stress  intensity  factors  when  Implemented  numerically,  and  either  can 
be  made  quite  computationally  efficient  when  used  together  with  a  super¬ 
position  procedure  or  regularizing  approach  such  as  {4] . 
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APPENDIX  ----- . 

Here  for  completeness  we  furnish  details  of  the  complementary  singu¬ 
lar  eigenfunctions  required  in  the  definitions  of  H_  and  H_I  in  (5),  (6). 
These  fields  are  most  readily  expressed  in  cylindrical  polar  coordinates 
r,  e  (Fig.  1).  For  the  Integral  we  then .have 

o*r  -  -kV^KX  +  3)cos(X  +  1)6  -  8  cos(X  -  1)6] 

a*  -  K*r“*-1[(X  -  l)cos(X  +1)6-8  cos(X  -  1)6] 

v“ 

e*  -  -K*r"X_1I(X  +  l)sin(X  +1)6-8  *in(X  -  1)6] 

To 

•Such  a  strategy  offers  computational  savings  with  finite  element  analysis 
of  about  a  factor  of  four  in  terms  of  operational  counts  and  around  a  fac¬ 
tor  of  two  in  storage  requirements.  See  [2],  Eqn.  (10)  et  seq.  for  specifics 
of  the  decomposition. 
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Ur  "  '2WX  +  *)cos(*  +  1)6-6  cos(X  -  1)6]  J 

*  -X 

Ue  "  I(X  "  K)*in<X  +  1)6-6  »inU  -  1)6] 

_*  m  /(8it)yX(X  sin2  a  +  sin2  la) 

(l  +  k)(s1t>  2a  +  2a  cos  2Xa) 

in  which 

l  .  !+  a  m  X2  -  1 

*  X-cos-2a""+"Tos’"Jxa 


and  wherein  p  la  the  shear  modulus,  and  k  “  3  -  4v  for  plane  strain  and 
(3  -  v)/(l  +  v)  for  plane  stress,  v  being  Poisson's  ratio.  For  the  H 
Integral  we  have  11 


o 


** 

rr 


**  —  X— l 

-K  r  a[B  sln(X  -  1)9  - 


(X  +  3)sln(X  +  1)9] 


°99  *  K*VX_1[8  sln(X  -  1)9  -(X  -  l)sln(X  +  1)9] 

°r9  -  K*VX_116  cos(X  -  1)9  -(X  +  l)coa(X  +  1)9] 
**  .x 

u**  -  [6  aln(X  -  1)9  -(X  +  Oeln(X  +  1)9] 

^  **  -x 

ue  -  -  [6  co.(X  -  1)9  -(X  -  k)cos(X  ♦  1)9] 

**  m  /(8n)yX(X  sin2  a  -  sin2  Xa) 

(1  +  «c)(sln  2a  -  2a  cos  2Xa) 

In  which  now 

X  .  .  _  X2  -  1 

*  X  cos  2a  -  cos  2Xa 


Like  fields  for  the  node  111  Integral  (-K^j)  can  be  found  In  [2]. 
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CORRIGENDUM 

Correction:  "A  Remark  on  the  Determination  of  Mode  I  and  Mode  II  Stress 
Intensity  Factors  for  Sharp  Re-Entrant  Corners, *'  G.B.  Sinclair,  Inter¬ 
national  Journal  of  Fracture  27  (1985)  R81-R85. 

Equations  (5)  and  (6)  on  pg  R82  should  read: 


Hi  “/(aij  “*  '  °*i  Ui)V8 

hh  -  °Vi  wb 


The  paragraph  following  Eqn.  (7)  on  pg  R82  should  read: 

Turning  to  their  application  in  mixed  mode  problems  which  have  both  sym¬ 
metric  and  antisymmetric  singular  fields  present,  we  first  focus  on  the 
use  of  H  to  determine  K  in  such  cases.  In  view  of  the  path  indepen¬ 
dence  ofTl  established  for  any  problem's  elastic  fields  in  [2],  the 
general  contour  I  can  be  exchanged  for  the  here  more  convenient  path  T 
(Fig.  1) - 

Equation  (8)  is  split  and  occurs  on  the  bottom  of  pg  R82  and  the  top  of 
pg  R83.  The  arguments  for  the  last  displacement  field  occurring  on  pg 
R83  should  read: 


UjUj.-Xj) 


s# 

r>.v' 


Ve  regret  any  inconvenience  created  for  our  readers  by  these  typographical 
errors  -  Ed. 
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Abstract 

A  number  of  current  practices  in  fracture  mechanics  which  use  quantities  near  a  crack  tip  to  make  conclusions 
about  response  at  the  crack  tip  itself  are  examined  Specifically  these  include:  stress  and  displacement  matching 
to  estimate  stress  intensity  factors,  monitoring  local  stress  and  strain  values  to  predict  fracture,  and  both  crack 
opening  angle  and  crack  opening  displacement  as  fracture  criteria  By  means  of  a  pair  of  counter  applications,  all 
of  these  procedures  are  demonstrated  to  have  the  potential  of  leading  to  completely  incorrect  conclusions.  An 
understanding  of  m-hat  causes  this  inadequate  performance  then  indicates  that  such  procedures  may  be  unreliable 
in  general  and  prompts  suggestions  as  to  alternatives. 


1.  Introduction 

In  fracture  mechanics  today  there  are  a  number  of  procedures  which  in  essence  draw  on 
field  quantities  in  the  vicinity  of  a  crack  tip  to  infer  what  is  happening  right  at  the  crack 
tip.  This  paper  considers  two  classes  of  such  procedures:  local  fitting  methods  for 
determining  stress  intensity  factors,  and  local  fracture  criteria  for  predicting  fracture.  The 
intent  is  to  examine  the  reliability  of  these  approaches. 

Local  fitting  methods  for  calculating  stress  intensity  factors  at  cracks  arose  out  of  a 
need  to  extract  this  parameter  from  numerical  analyses,  particularly  finite  element 
analyses.  One  of  the  first  discussions  as  to  how  best  to  undertake  such  exercises  is  that  of 
Chan  et  al.  [1],  though  certainly  the  approach  was  in  use  prior  to  [1]  even  if  somewhat 
informally.  Moreover  the  approach,  together  with  spin  offs  such  as  the  nodal  force 
technique  of  Raju  and  Newman  (2],  continues  to  be  used  today  in  conjunction  with  both 
finite  element  and  numerical  integral  equations  analyses,  as  is  evident  in  several  papers  in 
recent  conferences  (3,4].  All  of  these  methods  entail  matching  near  but  not  at  the  crack  tip 
to  estimate  the  stress  intensity  factor  there;  the  question  considered  here  is  how  reliably 
can  the  attendant  extrapolations  be  carried  out. 

Local  fracture  criteria  attempt  to  predict  fracture  by  checking  some  quantity  in  the 
vicinity  of  the  crack  tip.  Among  the  more  natural  choices  of  quantity  to  this  end  are 
measures  reflecting  stresses  and  strains  at  the  crack  tip.  Stresses  are  usually  used  in  elastic 
analyses  and  typically  in  complex  configurations  (e.g.  for  failure  in  composites  as  in 
Chamis  [3],  and  for  biomedical  applications  as  in  Valliappan  et  al.  [6]).  *  Strains  are 


*  There  do  exist,  though,  examples  of  the  direct  use  of  local  stresses  in  elastoplastic  treatments,  tee  for  instance 
MiUer  et  al.  (7) 
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normally  preferred  when  significant  plastic  flow  accompanies  fracture  (e  g  Newman  [8]. 
Belie  and  Reddy  [9],  Kjm  and  Hsu  [10]).  Other  quantities  employed  in  this  role  are  the 
crack  opening  displacement  of  Wells  [11]  and  the  crack  opening  angle  of  Andersson  [12], 
the  former  having  gained  sufficient  acceptance  to  merit  a  British  Standard  [13]  to  govern 
its  measurement.  *  All  of  these  criteria  involve  comparisons  made  away  from  the  crack  tip 
with  a  view  to  gauging  what  is  happening  there,  the  issue  of  concern  here  is  the  certainty 
with  which  such  comparisons  are  connected  to  fracture  at  the  crack  tip. 

We  begin  our  assessment  in  Section  2  by  specifying  the  local  fitting  methods  and  local 
fracture  criteria  for  testing/calibrating  and  thereafter  using  on  simple  elastic  problems  in 
Section  3.  The  conceptual  applications  in  Section  3  in  fact  have  known  solutions,  thereby 
furnishing  demonstrations  of  local  procedure  performance.  These  sample  evaluations  in 
turn  motivate  a  discussion  of  the  more  general  use  of  local  procedures  and  some 
suggestions  as  to  alternatives  in  Section  4.  The  paper  then  ends  with  some  concluding 
remarks. 

2.  Procedures  examined 

Here  we  first  define  a  stress  intensity  factor  and  describe  two  local  fitting  procedures 
aimed  at  estimating  it  -  one  using  stress,  the  other  displacement.  We  then  specify  some 
local  fracture  criteria  which  are  based  on  stress,  strain,  crack  opening  angle,  and  crack 
opening  displacement. 

To  fix  ideas  consider  a  cracked  elastic  plate  which  is  thick  in  the  out-of-plane  direction 
so  that  a  state  of  plane  strain  obtains  (Fig  1).  For  the  in-plane  directions,  we  let  (jr.  » )  be 
rectangular  cartesian  coordinates  having  origin  O  at  the  crack  tip  and  x-axis  aligned  with 
the  crack.  To  further  assist  the  development,  we  also  let  (r,  0)  be  cylindrical  polar 
coordinates,  sharing  the  same  origin,  and  related  to  the  rectangular  coordinates  by 

jr  “  r  cos  6,  y-r  sin  0.  (1) 

for  0 <  r<  x.  -v  <d<ir.  Next  we  assume  the  plate  to  be  under  symmetric  loading 
having  resultant  force  per  unit  thickness.  P,  acting  transverse  to  the  crack.  Then  the  only 
stress  intensity  factor  present  is  the  mode  1  factor  A,.  defined  by 

A |  «•  lim  \/2rrr OfU(Of)  on  #-0,  (2) 

r  —  0 

where  a,  is  the  elastic  normal  stress  in  the  ^-direction  and  V  is  the  unit  step  function,  here 
taken  as  being  one  for  positive  arguments  and  aero  otherwise. 

We  now  wish  to  draw  on  local  quantities  to  form  estimates  of  the  stress  intensity  factor. 
A,.  One  straightforward  stress  estimate  is  suggested  by  definition  (2)  itself  and  merely  sets 

*i-v2^b*#L  on  •-0.  (3) 

wherein  r0  is  some  “small'’  distance  from  the  crack  tip  and  the  bar  atop  o0  distinguishes  it 
as  being  whatever  value  is  found  by  the  numerical  solution  technique  adopted.  An 
analogous  displacement  estimate  is 

*/“TT on  6-rr.  (4) 

Here,  ji  is  the  shear  modulus,  a  Poisson's  ratio,  and  -ut  on  0- w  the  numerically-de- 

*  Crack  opening  displacement  mas  also  independently  introduced  by  Cottrell  |14)  to  the  effect  a  somewhat 
different  objective,  that  of  dassifyin|  brittle  versus  ductile  fracture  -  see  Burdekin  |1S)  for  a  recent  review  of  ns 
role  in  fracture  mechanics 
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Figure  1.  Mode  I  crack  configuration  and  coordinates. 


termined  opening  displacement  of  the  upper  crack  flank.  Again  in  (4).  for  simplicity,  we 
select  r  -  r0  as  our  “small"  distance,  a  policy  we  continue  throughout  the  remainder  of 
this  discussion  without  loss  of  generality.  The  estimates  in  (3),  (4)  both  fit  the  first  singular 
eigenfunction  for  a  symmetrically  excited  crack  to  the  entire  crack  fields  at  a  station 
supposedly  close  enough  to  ensure  that  these  singular  fields  dominate.  Thus  while  (3).  (4) 
do  not  represent  the  most  sophisticated  of  local  fitting  procedures  for  finding  K,.  they  do 
contain  the  essence  of  the  rationalization  for  all  such  approaches.  Consequently  they  are 
adequate  for  demonstration  purposes  here;  subsequently  we  review  other,  more  complex, 
local  fitting  methods. 

We  now  turn  to  fracture  criteria  based  on  local  quantities.  For  a  stress  criterion  we 
choose  perhaps  the  most  obvious,  namely  that  the  tensile  stress  ahead  of  the  crack,  near 
but  not  at  the  crack  tip.  attain  a  critical  value  for  fracture.  That  is. 

o#Ur0-Oc  on  *“0.  (5) 

for  fracture,  wherein  oc(  >  0)  is  the  critical  stress  determined  on  a  suitable  calibration 
problem.  An  analogous  strain  criterion  exists.  Instead,  with  a  view  to  examining  a  greater 
variety  of  local  fracture  criteria,  we  consider  the  more  involved  condition  of  Newman  |8). 
In  Newman  (8],  the  crack-tip  strain  measure  is  assembled  from  a  finite-element  analysis. 
The  result,  denoted  c,  here,  is  the  average  of  the  element  strains  in  the  >>-direetion. 
taken  over  all  the  elements  connected  to  the  node  at  the  crack  tip.  Accordingly  for  fracture 
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we  have  (Fig.  2) 

3 

I  (6) 

where  tc(  >  0)  is  a  calibrated  critical  strain  and  r,  -  r2  -  v^5  r0  /3.  r,  -  \/2  r0/3,  0,  -  cot'  '2. 
d2«tan'’2,  0} « 3w/4.  The  sample  finite  element  grid  of  Fig.  2  also  suffices  for 
computing  the  crack  opening  angle,  a.  of  Andersson  [12].  this  being  the  angle  subtended 
by  the  closest  node  to  the  crack  tip.  The  associated  fracture  criterion  has  that,  at  fracture. 

««ar.  a-2tan'*(  -w#U,/r0)  on  (7) 

with  <t,(>0)  being  a  calibrated  critical  angle.  A  similar  evaluation  of  the  crack  opening 
displacement.  6,  of  Wells  [11],  leads  to  the  fracture  criterion  that  has.  at  fracture. 

5«$r,  5-  -2u#|,_,0  on  0-ir.  (8) 

6r(  >  0)  being  a  calibrated  critical  displacement.  Actually  6  is  usually  evaluated  as  the 
crack  opening  displacement  at  r  -  rY  »  Af/hra}.  where  <r>  is  the  uniaxial  yield  stress: 
here  we  take  ey  to  be  such  that  this  station  too  coincides  with  r«r0.  Equation  (8) 
completes  our  specification  of  the  set  of  fracture  procedures  based  on  local  quantities  to 
be  appraised  next. 


3.  Test /calibration  problem:  applications 

In  this  section  we  start  by  formulating  a  test/calibration  problem  and  a  pair  of  problems 
to  serve  as  “applications”.  Then  we  exhibit  the  closed  form  solutions  to  all  three  problems, 
employ  the  First  of  these  to  test  the  local  fitting  methods  and  calibrate  the  local  fracture 
criteria,  and  thereafter  use  the  last  of  the  solutions  to  assess  performance. 

The  problems  to  be  treated  are  drawn  from  the  class  of  symmetric,  plane  strain,  crack 
problems  outlined  in  the  previous  section.  We  do  this  conveniently  by  confining  attention 
to  the  cracked  circular  plate  of  radius  R  (Fig.  3)  and  hence,  by  virtue  of  symmetry,  need 
only  to  consider  the  semicircular  region  3t,  defined  by 

{(r,d)|0<r</l,0<«  <»r}.  (9) 

Throughout  this  region  we  seek,  in  general,  the  stresses  or.  a,,  r*  and  displacements  ur,  u, 
as  suitably  smooth  functions  of  r,  9  meeting  the  following  requirements.  All  three 
problems  are  to  satisfy:  the  plane-strain  stress  equations  of  equilibrium  in  the  absence  of 
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Figure  3.  Cracked  plate  geometry  for  lest  problem  and  applications. 


body  forces,  namely 

re, ,  +  r.«  t  +  o.  -  o.  <■  0, 

(10) 

ot  i  +  r-T^  ,  +  2t^  —  0. 

on  #.  where  subscripts  preceeded  by  a  comma  denote  differentiation  with  respect  to  that 
variable:  the  plane-strain  stress-displacement  relations  for  a  linear  elastic,  homogeneous 
and  isotropic,  plate  with  shear  modulus  pi  and  Poisson's  ratio  »  «*  1/4.  via.. 

•r -M  [3«,., +  r_,(ii, +  «,.,)]. 

®«-M[3r',(M,+  «## )  +  «,,].  (11) 

on  the  stress-free  crack  flank  conditions  which  set 

on  9-v,  (12) 

for  0  <  r  <  R,  and  the  symmetry  conditions  ahead  of  the  crack  which  have 

m#“0,  r^»0  on  #«0,  (13) 

for  0  <r<R.  In  addition  each  problem  individually  must  comply  with  its  applied  stress 
conditions  on  the  circular  boundary.  These  set 


-  F  V 

®,  "  IB  2. 

P  r' 

Tf LB  *-> 


a,  cos—  on  r—R, 

«  n9  n 

e,,sin-y  on  r-/?. 


for  0  <  9  <  ir(  P  >  0),  where,  for  the  test  problem. 
*■8,  M|  -  5,  O]  -  - 1 .  h,-b3-l. 


are  the  only  a,  #  0,  bH  *  0,  and.  for  the  two  applications,  all  nontrivial  constants  are  given 
in  Table  1.  wherein  p  -  R/r0(  >1). 

While  not  strictly  needed  in  the  mathematical  formulation  of  the  foregoing  type  of 
elasticity  problem,  we  adjoin  the  following  regularity  requirements  on  the  displacements  in 
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Tahlf  1.  Boundary  condition  consume  for  the  applications 


Constant 

Application  t 

Application  2 

A 

8fl,,(p  -2+j2 

8<p3-l>: 

-7.374.036p171 

-  5  969.792  p1 71 

0 

5 

a  2  b: 

0 

15  690.91 8p'7j 

u  y 

10(9  -4&  >p!/<9 

-  1  -  10pJ 

Uk 

2(10+  i/2  lp'/7 

0 

u*  -  K 

0 

47.072.75 4p’  : 

«■* 

3<flj-  P,,/7) 

-  3p:  ( 1 0  +  ,9p: ) 

-  <>» 

0 

10.008.616p’  : 

<1 4 

-  $fl, 

0 

fln-  — 

*  1 1  *  6/2  lp‘/7 

6.V* 
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all  three  problems  to  emphasize  their  physical  admissibility:  that  the  displacements  at  the 
crack  tip  be  bounded,  that  is 

u,«0(l).  «#-0(l)  as  r-»  0.  (16) 

on  £.  where  O  is  the  large  order  symbol:  and  that  the  crack  opening  displacement  be 
nowhere  negative,  i.e.. 

-«/#>0  on  6~r.  (17) 

for  0  <  r  <  R. 

We  postpone  for  the  moment  a  discussion  of  how  the  test  problem  and  two  applica¬ 
tions  were  arrived  at  together  with  their  solutions  and  merely  exhibit  the  latter  at  this 
point.  For  all  three  problems  the  complete  solutions  admit  lo  representation  by  finite 
sums  of  eigenfunctions  for  the  crack,  viz.. 
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on  &.  For  each  problem  respectively,  the  solutions  have  k  as  previously  ((IS).  Table  1) 
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Table  2.  Coefficients  in  the  solutions  for  the  applications 


Coefficient 

Application  1 

Application  2 

fi 

0 

1 

'u  2(: 

-1.843.S09 

-1.492.448 

3  ; 

'U  <4 

0 

1.845.459 

10(9-4^1/49 

-10 

'll  ■‘V 

0 

1.251,077 

X 

fl»<T 

-2(10+  \/2  1/7 

0 

ft»‘  * 

901  +  6^1/49 
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while  the  only  nontrivial  coefficients  in  the  above  are:  for  the  test  problem, 

c.-l:  (19) 

and  for  the  two  applications,  those  c„  #  0  in  Table  2. 

The  solutions  in  (18).  (19).  Table  2  may  be  verified  directly.  Substituting  the  forms  in 
(18)  into  (10).  (11)  establishes  that  the  field  equations  are  complied  with,  and  inspection  of 
(18)  shows  that  the  crack  flank  and  symmetry  conditions  (12),  (13)  are  met.  Evaluating  o,, 
of  (18)  on  r  «  R,  using  (19),  Table  2,  further  shows  that  the  boundary  conditions  (14), 
(15).  Table  1  are  fulfilled.  Finally  ,  checking  the  displacements  in  (18)  under  (19)  or  Table 
2  reveals  that  the  regularity  requirements  (16),  (17)  are  satisfied. 

In  applying  the  local  fracture  procedures  to  the  three  problems  we  begin  by  assuming 
that  their  numerical  treatments  have  been  refined  to  the  extent  that  the  exact  answers  are 
recovered  to  all  intensive  purposes,  i.e.  8,  -  ot ,  etc.  This  is  the  ultimate  of  situations  from 
an  analy  sis  point  of  view  and  allows  us  to  focus  on  whatever  errors  are  introduced  by  the 
procedures  alone. 

Invoking  this  assumption  we  next  test  our  local  fitting  methods  before  applying  them,  a 
precaution  commonly  undertaken  in  practice.  Using  the  local  fits  for  the  stress  intensity 
factor  in  (3),  (4)  with  p  «  1/4  on  (18)  when  (19)  holds  yields,  for  both  methods. 


For  comparison,  taking  the  defining  limit  of  (2)  in  (18),  (19)  gives 


*/- 


4P  [2^ 
k\R 


(21) 


Hence  the  local  fitting  values  of  AT,  are  exact  and  the  performance  of  the  two  methods  (3), 
(4)  on  the  test  problem  is  perfect. 

In  order  to  calibrate  our  local  fracture  criteria  we  need  to  determine  the  load  to  fracture 
in  our  standardizing  problem.  We  do  this  by  regarding  fracture  as  occurring  when  the 
stress  intensity  factor  takes  on  a  limiting  value  of  Kr%  the  fracture  toughness  for  the  plate. 
It  follows  from  (21 )  that  the  critical  load  in  the  calibration  problem,  Pe,  is 


Thus  increasing  P  to  P(  and  evaluting  the  stress  in  (5),  the  strain  measure  of  (6)  *,  the 


*  For  this  derivation  we  note  that  straightforward  manipulation  leads  to 

«,  "  T~7~rs  I  r.r..-i./j[(2-._(_)'2)«*(,-2)f^(«-J)«»(n— 6)|] 
ipM'R  .-u. ...  1  ‘  ^ 

as  the  norma)  strain  in  the  direction  to  be  calculated  at  the  set  of  points  in  (6)  et  seq..  and  thereafter  averaged 
to  produce  . 
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crack  opening  angle  (COA)  of  (7).  and  the  crack  opening  displacement  (COD)  of  (8). 
using  (18).  (19)  gives,  respectively. 

=  1  _  *Pr  K, 


k\r0R  y'2  vr0 
4  Pck'  K,k' 


on  6  -  0, 


kn\r0R  fik2wr0 
l )P ", 


a  -  2  tan-1 


’2  tan" 


3  K, 


(23) 


pk^r0R  2^2vr0 

12  Pc  [r^  2Kr  fr^ 


5‘7 ft/s-irVs-8- 

wherein  p  has  been  set  to  1/4  and  k'  —  0.303,299. 

With  our  local  fitting  methods  tested  and  local  fracture  criterion  calibrated  we  can 
apply  them  to  the  two  applications.  Introducing  (18),  together  with  the  coefficients  of 
Table  2,  into  (3).  (4).  (5),  (6).  (7)  and  (8)  with  ►  -  1/4  and  when  F-  P(  of  (22).  with  k 
therein  now  taking  on  the  respective  values  in  Table  1,  furnishes  the  results  displayed  in 
Table  3.  As  opposed  to  normal  engineering  practice,  we  can  also  determine  the  corre¬ 
sponding  exact  answers  for  the  stress  intensity  factors  in  our  “applications"  from  (2).  (18) 
and  Table  2.  and  thereby  check  the  values  of  when  P  -  Pc.  These  last  results  are  also 
included  in  Table  3. 


Table  3.  Comparison  of  procedure  outcomes  with  actual  results 


Activity 

Application  ] 

Application  2 

A  ,  detcnniiution 

Suns  Tit 

*•-¥/? 

o 

1 

Displacement  fit 

*1 

1 

o 

Exact  answer 

0 

Fracture  prediction 

Local  stress 

* 

• 

■ 

J* 

* 

i 

o' 

1 

O 

Local  strain 

iy  “  €C 

f,-0 

COA 

a-0 

COD 

1-0 

Exact  answer 

F-A-*t- 0 

Several  comments  on  the  results  in  Table  3  are  in  order.  For  the  first  application,  the 
stress  intensity  factor  in  reality  is  zero  irrespective  of  load  level  P  so  that  the  stress  right  at 
the  crack  tip  is  also  always  zero  and  fracture  there  is  not  a  possibility  (at  least  not  prior  to 
fracture  elsewhere):  the  local  fitting  procedures  on  the  other  hand  estimate  K ,  as  being 
dependent  on  P  and  consequently,  when  this  load  is  sufficiently  large  (P  -  Pc  in  fact),  all 
the  local  fracture  criteria  predict  fracture  at  the  crack  tip.  Hence  the  procedures  consid¬ 
ered  give  rise  to  conservative  but  gross  errors  in  Application  1.  For  the  second  application, 
the  stress  intensity  factor  is  not  zero  and  increases  linearly  with  load  so  that  there  is 
indeed  a  load  level  at  which  fracture  occurs:  the  local  fitting  procedures  in  contrast  give 
no  K  |  whatsover,  independent  of  the  value  of  P,  and  companion  predictions  of  no 
fracture  for  any  load.  Hence  the  procedures  give  rise  to  nonconservative  gross  errors  in 
Application  2.  In  all  the  results  in  Table  3  represent  a  dear  demonstration  of  the  potential 
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unrealiability  of  (he  fracture  procedures  of  Section  2  -  we  look  for  means  of  controlling 
such  unreliability  next. 

4.  Some  alternatives 

Here  we  start  by  examining  the  nature  of  the  breakdown  of  the  two  local  fitting  methods 
used  to  estimate  Kt.  In  the  light  of  this  appraisal  we  review  other  local  fitting  methods  and 
advance  suggestions  for  improved  methods  of  extracting  stress  intensity  factors.  Finally 
we  discuss  what  remedial  action  can  be  taken  for  the  local  fracture  criteria. 

The  reasons  underlying  the  failure  of  the  X,  estimation  via  local  fitting  in  the  preceding 
section  are  as  follows.  First,  the  fit  must  be  made  away  from  the  crack  tip  with  r0  *  0, 
since  the  stresses  can  be  singular  and  thus  not  finable  at  r0  «  0  while  the  displacements 
are  zero  there  leaving  nothing  to  fit.  Second  and  as  a  consequence,  other  eigenfunctions 
for  the  crack  that  are  not  being  fitted  can  participate  causing  incon-ect  answers.  Indeed 
this  is  the  effect  that  the  applications  in  Section  3  were  constructed  to  produce.  That  is. 
symmetric  eigenfunctions  for  the  crack  other  than  the  singular  one  being  fitted  were  added 
so  as  to  generate  the  erroneous  K,  estimates  of  Table  3,  then  whatever  stresses  the 
combination  so  devised  realized  on  the  circular  boundary  were  taken  as  the  prescribed 
stresses  there.  • 

At  first  glance  it  might  appear  a  simple  matter  to  design  a  local  fitting  method  which 
does  not  succomb  to  such  contrivances.  Perhaps  the  simplest  candidate  stategy  to  this  end 
is  to  move  the  point  fitted  closer  to  the  crack  tip;  such  an  approach  however  fails  if  r0  in 
the  counter  applications  is  then  adjusted  to  coincide  with  the  new  point  (which  it  can  be). 
A  more  sophisticated  possibility  is  the  procedure  in  which  Pietro,  on  8  -  0  is  plotted  as  a 
function  of  r,  a  straight  line  fitted  to  the  data,  then  the  line  extrapolated  to  r  «  0  to 
estimate  X,.  This  type  of  extrapolation  technique  in  effect  matches  the  first  and  third 
eigenfunctions  for  the  crack  and  would  appear  to  enjoy  some  popularity  in  present  day 
practice.  There  are  a  number  of  ways  in  which  the  details  of  this  and  like  extrapolation 
methods  may  be  implemented  with  varying  results,  but  typically  such  an  approach  does 
meet  with  more  success  on  the  given  applications  than  that  reported  in  Table  3.  In  fact,  on 
knowing  the  answers,  it  is  possible  to  adjust  the  technique  to  recover  the  correct  results 
exactly.  For  such  a  fine  tuned  method  though,  it  is  then  possible  to  devise  a  further  new 
application  on  which  it  does  not  work  by  superimposing  even  more  eigenfunctions.  By 
way  of  illustration  of  this  son  of  more  extensive  counter  application,  consider  Fig.  4.  **  In 
Fig.  4,  /2 vrot/K  on  8  -  0  is  given  as  a  function  of  r/R  on  a  set  of  17  points;  here  K  is  a 
nondimensionalizing,  but  otherwise  arbitrary,  constant.  A  least  squares,  straight  line 
through  these  data  points  passes  through  the  origin,  yielding  K\  -  0.  In  actuality.  A, «  X. 
illustrating  the  unrealiability  of  this  approach.  Further,  gathering  more  data  on  1/9  <  r/R 
<  1  and  matching  does  not  improve  the  incorrect  estimate  of  Xj.  Furthermore,  similarly 
erroneous  results  can  occur  when  using  the  displacement  counterpart  -  indeed,  there  is  a 
counter  application  which  has  exactly  the  same  data  as  that  depicted  in  Fig.  4  but  with 
"•Jlvrot/K  on  8 -  0”  there  exchanged  for  ijn/lr  fiut/K(  1  -  r)  on  8-v'\  thus 
producing  Jf,  ■  0  despite  the  fact  that  X,  —  X. f  And  this  is  the  situation  that  prevails  in 
general.  That  is,  irrespective  of  how  apparently  refined  a  local  fitting  method  is  and  what 
it  fits,  once  its  specifics  have  been  decided  on  it  is  then  always  possible  to  generate  an 
example  on  which  it  proves  to  be  comprehensively  inadequate.  The  reason  for  this  is  that 


*  la  this  protest  the  choice  of  »  - 1/4  it  amply  made  to  ((dilute  the  construction  and  not  of  any  afnificance 
m  itself 

••  See  Sinclair  {16]  for  details  of  this  example. 

1  See  {16]  for  details  of  this  second  example  and  others. 


12 


O  B  Sinclair 


f/R 

Figure  4  Product  of  ttreu  ahead  of  the  crack  times  the  square- root  of  distance  from  the  crack  in  a  more 
extensive  counter  application. 

any  local  fitting  procedure  must  match  a  finite  number  of  quantities  on  a  finite  set  of 
points  and  thereby  can  only  fit  a  finite  number  of  eigenfunctions;  thus  any  local  fitting 
procedure  leaves  an  infinity  of  eigenfunctions  unfitted  whose  participation  factors  can  be 
adjusted  to  ensure  its  downfall. 

A  key  question  that  arises  at  this  point  therefore  is  who  has  the  final  say.  the  analyst 
making  adjustments  to  his  local  fitting  method  so  that  it  does  produce  good  results  or  the 
rascal  generating  the  problems  on  which  the  given  method  does  not  work.  Unfortunately 
in  practice  the  latter  always  has  the  option  of  going  last.  This  is  because  in  practice  we 
must  ultimately  stop  checking  our  procedure  on  test  problems  and  apply  it.  Then,  since 
applications  are  problems  for  which  we  do  not  know  the  answers,  we  cannot  tell  with 
certainty  how  well  our  procedure  really  works,  and  as  a  result  we  have  to  tolerate  the 
possibility  of  it  giving  erroneous  estimates,  maybe  even  grossly  erroneous  estimates. 

Logically  one  only  needs  the  possibility  of  such  results  in  an  application  to  establish 
that  a  method  is  unreliable.  In  practice,  howeveT,  one  can  perhaps  rationalize  still  using 
the  method  if  one  can  nssign  a  sufficiently  low  probability  to  the  occurrence  of  the 
disastrous  example  •  the  counter  application  being  pathological  to  a  degree.  Certainly  the 
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counter  applications  provided  in  Section  3  do  not  appeal  as  physical  problems,  being 
constructed  more  towards  dramatically  underscoring  how  much  in  etror  the  selected  local 
fitting  methods  and  fracture  criteria  could  be  Even  so,  the  two  particular  counter 
applications  are  one  of  an  infinity  of  pairs  that  could  be  devised  to  ensure  that  the 
methods  in  (3).  (4)  give  the  same  completely  incorrect  estimates  as  in  Table  3  Moreover 
there  are  other  such  sets  of  pairs  that  give  rise  to  estimates  that  are  more  subtly  in  error 
than  those  of  Table  3  but  are  nonetheless  significant!)  inaccurate  And  such  is  the  case  for 
any  local  fitting  method,  no  matter  how  apparently  refined  -  one  can  always  construct  an 
indefinitely  large  number  of  applications  on  which  it  fails.  It  therefore  appears  unlikely 
that  the  probability  of  encountering  a  problem  application  in  engineenng  can  safely  be 
taken  as  being  negligible. 

More  specifically  in  this  regard,  consider  the  following  limited  simulation  of  perfor¬ 
mance  for  the  simple  stress-fitting  method  of  (3)  on  two  problems  that  could  reasonably 
be  viewed  as  being  of  practical  consequence.  In  accordance  with  normal  practice,  we  begin 
our  simulation  by  testing  our  method.  As  a  test  problem  we  take  the  classical  Griffith 
geometry  of  a  center-cracked  infinite  plate  and  consider  the  case  when  the  crack  is  opened 
by  a  uniform  far-field  tension.  Then  we  adjust  r0  so  as  to  ensure  acceptable  results  which 
here  we  take  as  being  up  to  5$  in  error.  Picking  r0  -  a/ 14.  a  being  the  semi-crack  length, 
as  a  moderately  close  point  at  which  to  obtain  accurate  stress  values  by  either  numerical 
or  experimental  methods,  we  find  on  assuming  negligible  errors  in  the  stress  determination 
that 

£(-1.03*1.  (24) 

wherein  A,  is  the  exact  value.  Equation  (24)  represents  a  satisfactory  and  encouraging 
result.  Now  as  an  “application”  we  take  the  Griffith  crack  opened  by  a  uniform  pressure. 
Applying  (3)  with  r0  -  a/14  as  before  then  yields 

£,-0.67*;.  (25) 

Equation  (25)  represents  an  unsatisfactory  nonconservative  result,  showing  that  the 
probability  of  meeting  a  problem  application  in  practice  for  this  simple  method  is  not 
zero.  • 

Of  course,  as  remarked  earlier,  it  is  possible  to  adapt  ones  local  fitting  method  to 
overcome  the  difficulty  experienced  in  this  last  elementary  example.  However,  while 
increasing  ones  efforts  considerably  when  applying  local  fitting  procedures  by  performing 
the  matching  in  a  variety  of  ways  at  a  number  of  locations  does  tend  to  increase  the 
chances  of  detecting  spurious  results,  it  does  not  appear  reasonable  at  this  time  to  assume 
that  any  such  combination  of  techniques  reduces  the  probability  of  unacceptably  incorrect 
results  to  a  negligible  level  in  all  situations  of  practical  interest.  Continuing  to  use  local 
fitting  methods  to  determine  stress  intensity  factors  in  engineering  would  thus  seem  to  be 
unjustified  and  unwise.  This  is  especially  so  since  there  are  available  quite  different 
approaches  which  are  free  from  the  potential  of  furnishing  useless  results  that  local  fitting 
methods  have,  and  which  can  be  implemented  with  no  more  effort. 

One  of  the  best  of  such  alternative  approaches  employs  path  independent  integrals, 
since  the  integral  operators  involved  are  essentially  orthogonal  to  all  of  the  eigenfunctions 
for  the  crack  except  the  singular  one.  That  is,  contributions  to  these  integrals  from  the 
nonsingular  eigenfunctions  can  be  shown  to  be  zero  analytically  **,  so  that  these  terms  are 

*  See  Efus  et  el.  (17)  for  further  discussion  of  the  effecu  of  the  presence  of  stresses  which  ere  independent  of  r, 
like  the  additional  hydrostatic  pressure  here.  See  also  Sinclair  and  Multan  [II]  for  a  further  example,  concerns nj 
the  standard  configuration  of  a  single  edge  notch  under  tension,  wherein  local  fitting  methods  have  been 
responsible  for  less  than  satisfactory  results 

••  See  e  g.  |19]  pp  1002.  1003  for  a  proof. 
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in  essence  eliminated  rather  than  merely  hoped- to-be-small  as  they  are  to  varying  degrees 
in  local  fitting  methods.  Two  such  integrals  are  the  J  integral  of  Eshelby  [20],  Rice  [21] 
and  Sanders'  integral  [22];  both  of  these  enjoy  physical  interpretation  as  the  energy  release 
rate  at  the  crack  tip  in  the  elastic  instance.  Other  such  integrals  are  those  of  Carpenter  [23] 
and  of  Sinclair  et  al.  [19];  these  last  lack  the  appeal  of  a  direct  physical  interpretation  but 
can  be  expected  to  perform  somewhat  better  numerically  since  they  involve  the  calculation 
of  single  terms  rather  than  products.  *  In  addition  there  are  also  methods,  like  that 
developed  by  Parks  [24],  which  effectively  employ  a  path  independent  integral  and 
consequently  realize  similar  performance. 

Turning  to  fracture  criteria  based  on  local  stresses  and  strains,  we  find  the  situation  is 
not  improved  in  terms  of  their  likelihood  of  furnishing  correct  information  Again  one  can 
conceive  of  more  complex  criteria  than  the  simple  stress  criterion  and  the  average  strain 
criterion  specifically  examined  here;  for  example,  the  condition  of  Belie  and  Reddy  [9] 
which  checks  the  maximum  strain  in  all  the  elements  near  the  crack  tip  to  decide  if 
fracture  occurs  there.  Again  too  one  can  construct  counter  applications  on  which  such 
approaches  fail;  for  example.  Application  2  actually  has  t,  «0  at  each  centroid  of  the 
elements  sharing  the  crack-tip  node  (Fig.  2).  so  that  Belie  and  Reddy's  criterion  applied  to 
these  elements  would  predict  no  fracture  despite  the  fact  that  crack-tip  fracture  occurs  in 
this  application  when  the  load  is  sufficiently  large.  Simply  put.  these  methods  cannot 
reliably  infer  what  the  physical  stresses  and  strains  are  near  the  crack  tip  given  an 
underlying  theory  which  can  lead  to  infinite  stresses  and  strains  at  the  crack  tip  in 
response  to  infinitesimal  loads.  Good  alternatives  here  await  the  development  of  stress 
and  strain  fields  for  cracks  whose  physical  relevance  is  unquestionable.  In  the  meantime, 
the  mathematically  singular  stress  and  strain  fields,  while  not  a  physical  reality,  are  a 
reality  of  analysis  which  must  be  faced,  even  in  the  elasto-plastic  instance.  **  At  this  point 
the  accepted  way  of  doing  this,  at  least  in  the  elastic  situation,  is  by  means  of  the 
associated  energy  release  rate,  whence  the  stress  intensity  factor.  Accordingly  local 
stress/strain  fracture  criteria  must  be  interpreted  as  trying  to  estimate  stress  intensity 
factors,  whether  wittingly  or  otherwise.  It  follows  that  such  practices  are  every  bit  as 
subject  to  the  uncertainties  of  explicit  estimation  via  local  fitting  methods  discussed 
previously,  and  we  recommend  that  alternatives,  such  as  those  suggested  earlier,  be  used  in 
their  stead. 

The  last  two  local  fracture  criteria  could  both  bi  regarded  as  using  measures  of  the 
loaded  crack  profile  and  thereby  avoiding  the  singular  stresses  and  strains  to  some  extent. 
Nonetheless  their  performance  in  the  elastic  regime  is  impaired  by  the  same  shortcoming 
that  all  of  the  preceding  local  procedures  have,  that  of  not  being  able  to  reliably  assign 
what  proportion  of  their  totals  are  due  to  benign  nonsingular  eigenfunctions  and  accord¬ 
ingly  what  part  is  not.  As  a  result  they  cannot  reliably  estimate  the  accepted  governing 
parameter,  the  stress  intensity  factor. 1  Furthermore,  being  fits  of  a  single  parameter  at  a 
single  point  in  effect,  they  cannot  even  be  readily  supplemented  to  reduce  the  chances  of 
extraneous  results.  In  the  light  of  this  inadequacy  of  both  measures  in  treating  elastic 
response  it  would  not  appear  reasonable  to  entertain  with  any  great  confidence  their  use 
to  treat  fracture  when  significant  plastic  flow  is  present,  the  end  to  which  they  were 
originally  proposed.  Especially  since,  physically,  elastic  response  must  precede  plastic 
irrespective  of  how  much  plastic  flow  ultimately  accumulates.  In  the  absence  of  a  sure 

*  Tim  has  in  reality  been  found  lo  be  the  cate  in  limited  numerical  experiments  undertaken  to  dale,  refer  118.19] 

**  See.  for  example.  Hutchinson  (25)  for  an  analysis  demonstrating  the  persistence  of  singular  behavior  within 
the  deformation  theory  of  plasticity. 

’  In  passing  we  observe  that,  in  view  of  the  uncertainty  in  the  elastic  situation  of  COD'S  relationship  to  A,, 
hence  to  J.  attempts  to  find  a  truly  general  simple  relation  between  i  and  J  in  elasto-plastic  instances  would 
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connection  to  a  physically  reasoned  explanation  of  fracture  then,  these  measures  remain 
inherently  unreliable  and  others  are  to  be  preferred.  The  best  of  these  at  this  time  would 
seem  to  be  the  stress  intensity  factor  itself  and  the  J  integral,  both  of  which  at  least  have 
the  attribute  of  being  clearly  related  to  the  energy  release  rate  in  the  special  case  of  pure 
elastic  response. 

5.  Concluding  remarks 

Fracture  mechanics  today  continues  to  be  faced  with  the  task  of  making  physically 
sensible  interpretations  of  the  nonphysical  singular  fields  w  hich  exist  at  crack  tips,  be  the 
situation  purely  elastic  or  elasto-plastic.  A  prerequisite  to  success  therefore  is  an  accurate 
assessment  of  the  degree  to  which  these  singular  fields  are  present.  The  procedures 
considered  here  -  local  fitting  methods  to  estimate  stress  intensity  factors  and  local 
fracture  criteria  to  predict  fracture  -  fail  to  do  this  in  a  way  which  can  be  shown  to  be 
reliable. 

The  basic  reasons  for  the  potential  unreliability  of  the  local  procedures  are  as  follows. 
First,  all  the  local  procedures  must  consider  quantities  near  but  not  at  the  crack  tip. 
.  Second,  at  such  stations  fields  other  than  the  key  singular  ones  can  contribute.  Third,  the 
extent  of  such  participation  cannot  be  either  completely  controlled  or  fully  accounted  for. 
As  a  consequence,  for  any  given  local  procedure  there  exist  problems  on  which  it  produces 
unacceptably  erroneous  results. 

The  question  then  arises  as  to  how  likely  is  one  to  come  upon  such  problem  problems 
in  practice:  the  answer  unfortunately  is  not  obvious.  For  local  procedures  based  on  a 
single  parameter  such  as  one  point  stress  or  displacement  matching  for  /fj.  a  critical  stress 
at  a  single  station,  COA,  and  COD.  the  likelihood  would  appear  to  be  quite  high.  For 
local  procedures  employing  a  multiplicity  of  quantities  such  as  some  of  the  more  extensive 
matching  techniques  to  estimate  stress  intensity  factors,  the  chances  would  seem  to  be 
reduced.  However,  for  the  best  of  these  approaches  in  this  regard,  the  probability  of 
invalid  results  cannot  be  shown  to  be  zero.  Consequently  in  real  applications  wherein  the 
answers  are  not  known,  one  cannot  be  sure  that  such  procedures  do  in  fact  find  them. 
That  is  not  to  say  that  local  procedures  may  not  work  on  occasions,  and  even  quite  well;  it 
is  just  to  say  we  cannot  be  certain  they  do.  Accordingly  local  procedures  must  be  viewed 
as  being  unreliable  at  this  time. 
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Resume 

On  examine  diverses  pratiques  courantes  en  Mecanique  de  Rupture  qui  font  usage  de  grandeurs  parametnques 
au  voisinage  de  l'exiremne  d  une  fissure  pour  nrer  des  conclusions  sur  ce  qui  se  passe  a  cette  extremne  merne 
Ces  grandeurs  torn  notamment:  les  valeurs  de  la  contrainte  et  du  deplacement  correspondant  a  un  facteur 
dimensne  de  contrainte  esttme.  l'enregistrement  des  contraintes  et  deplaccments  locaux  en  vue  de  predtre  la 
rupture,  et  l'angle  et  le  deplacement  d’ouverture  d’une  fissure  en  tant  que  criteres  de  rupture 
En  consider  ant  une  patre  d'appbcations  divergentes.  on  demontre  que  toutes  ces  procedures  nsquent 
d'entrainer  des  conclusions  compieiemem  incorrectes 

Une  analyse  des  causes  de  ces  madequanons  montre  que  ces  procedures  peuvent  se  reveler  peu  fiables  en 
general,  et  conduit  4  des  suggestions  de  procedures  alternatives. 
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FURTHER  EXAMPLES  OF  THE 
UNRELIABILITY  OF  STRESS  INTENSITY 
ESTIMATES  FOUND  VIA  LOCAL  FITTING 


INTRODUCTION 

The  inherent  unreliability  of  local  fitting  methods  for  determining  stress  intensity 
factors  is  established  in  [1,  23.  In  [1,  2],  examples  of  this  shortcoming  are  limited 
to  very  simple  fitting  procedures,  although  the  general  nature  of  source  of  the 
difficulty  is  discussed:  the  intent  in  this  brief  supplementary  report  is  to  furnish  the 
details  of  some  more  extensive  demonstrations.  We  begin  with  a  precise  statement 
of  four  chosen  problems,  then  exhibit  their  complete  closed-form  solutions*.  Using 
these  forms  we  next  extract  the  exact  values  of  the  corresponding  stress  intensity 
factors  and  compare  the  results  with  those  obtained  by  fitting  stresses  and 
displacements  using  several  eigenfunctions.  The  report  closes  with  a  brief  discussion 
of  some  of  the  implications  of  the  examples  regarding  even  more  involved  local 
fitting  procedures. 

FORMULATION  OF  PROBLEMS 

The  plane  region  of  concern  in  all  four  problems  is  the  cracked  circular  disk  of 
radius  R  «  9  (Fig.  1).  To  describe  this  geometry,  we  take  rectangular  cartesian 
coordinates  (x,  y)  with  origin  at  the  crack-tip  and  x-axis  aligned  with  the  crack  face, 
together  with  cylindrical  polar  coordinates  (r,  6)  related  to  (x,  y)  by 

x  •  r  cos  6.  y  «  r  sin  6  (0  S  r  <  oo,  -  V  <  S  S  V).  (1) 

Since  excitation  is  restricted  to  that  which  is  symmetric  about  the  x-axis,  we  can 
confine  attention  to  the  upper  half  of  the  disk,  R,  where 

R  »  {(r.  6)  |  0  S  r  <  R,  0  £  6  <  17}.  (2) 

We  seek  then,  the  stresses  ex,  o  g.  r(g.  and  displacements  uf,  u^,  as  sufficiently 
smooth  functions  of  r,  6  throughout  R,  satisfying  the  following  requirements.  The 
stress  equations  of  equilibrium  in  the  absence  of  body  forces. 


*tn  the  interests  of  brevity,  no  information  k  to  how  these  particular  problem*  were  constructed  is  given  in 
this  regard,  see  p.  2]  for  the  basic  approach. 


ait; 


3 


ffr  +  1  8rrg 
dr  r  dff 


Bffe,  Btt8  4 

dfl  dr 


(3) 


on  R.  The  plain  strain  stress-displacement  relations  for  a  homogeneous  and  isotropic, 
linear,  elastic  material  with  shear  modulus  y  and  Poission's  ratio  V, 
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on  R.  The  stress-free  crack-face  conditions  and  the  symmetry  conditions  ahead  of  the 
crack, 

°8  “  rx8  "  0  when  ^  *  W- 

rr#  *  0  when  0  *  (5) 

for  0<r<R.  The  stress  boundary  conditions  on  the  outer  circular  boundary  which 
prescribe 

K  ^  /<2n-1)0\ 

when ' '  "■ 

K  JL*  /(2n-1)0\ 

r  *  •  — wbsss,  >  b  sin\ - /  when  r  ■  R,  (6) 
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for  0<8<1T.  where  K  is  a  constant  (K>0)  and  a  ,  b  (n  »  0,  1....8)  are  given  in  Table  1 

n  n 

for  the  respective  problems  (1-4).  By  way  of  illustration,  Figs.2  display  sketches  of 
these  stress  distributions  for  Problem  1  (Figs.  2a  and  2b  share  the  same  scale).  And 
finally  the  regularity  requirements  which  insist  that  the  displacements  be  integrable  at 
the  crack  tip  and  that  there  is  no  overlapping  of  the  crack  flanks. 
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u.  *  0(1),  u^,  *  0(1)  as  r-»0, 
for  0^0^77 .  and 


£  0  when  8  *  77, 


for  0<r£R. 


TABLE  1:  BOUNDARY  CONDITION  CONSTANTS 


Constant 


Problem  1 


Problem  2 


Problem  3 


Problem  4 
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PROBLEM  SOLUTIONS 

For  all  four  problems  the  complete  solutions  can  be  expressed  as  finite  sums  of 
eigenfunctions  for  the  crack.  Explicitly  we  have 

K  4-  c  rn  5  1  1  3 


r  «  r—\  -21—  [( —  n)  cos(n  -  -)6  *  (n  — )  cos(n  ♦  -81. 
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a  >«T  [<n  ♦  -)  cos(n  -  -)8  -  (n  — )  cos(n  ♦  -)81, 
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where  c,  c  (n=0,1 . 6)  are  constants  having  values  for  the  respective  problems  as 

n 

given  in  Table  2.  That  (9)  and  Table  2  constitute  the  solutions  to  the  four  problems 
posed  earlier  can  be  verified  directly .  That  is,  the  forms  in  (9)  for  each  n 
individually  can  be  shown  to  satisfy  the  governing  field  equations  (3),  <4)  by 
substitution  and  the  crack-face  and  symmetry  conditions  (5)  by  inspection,  while  the 
combination  realized  on  the  boundary  r*R*9  using  the  values  in  Table  2  can.  after 
come  manipualtion,  be  shown  to  comply  with  the  prescribed  stresses  there  (6),  Table 
1.  In  addition,  inspection  of  (9)  shows  that  the  displacements  are  bounded  ar  r«0,  i.e. 
(7)  is  satisfied,  and  combining  the  forms  in  (9)  using  Table  2  gives  positive  crack¬ 
opening  displacements,  viz.,  (8)  is  met. 

TABLE  2:  SOLUTION  CONSTANTS 


Constant 


Problem  1 


Problem  2 


Problem  3 


Problem  4 


595,350 


c  -<l> 


0 

0 

651,657 

462.400 

121,615 


:•] 
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K,  ■  "" 

1  c 

Thus  from  Table  2  we  obtain  the  exact  values  of  K(  recorded  in  Table  3  for  each  of 
the  problems. 

As  out  local  stress-fitting  method  we  match  the  transverse  stress,  o g,  on  the  line 
of  symmetry  ahead  of  the  crack  (0=0)  at  the  five  radial  stations  r*1.  3,  4,  7,  9.  To 
do  this  we  take  the  five  eigenfunctions  for  the  crack  given  by  (9)  with  n=0,  1,  3/2,  2, 
5/2  therein**.  Then,  in  view  of  definition  (10).  we  must  solve  the  5x5  system  of 
equations 


k,  .  d,, .  <y3,!  *  i/  *  a/'*  *5^*91 


r  *  1.  3.  5,  7,  9, 

where  Kj  is  our  stress-fitting  estimate  of  the  stress  intensity  factor  and  dR  (n=1-4)  are 
estimated  participation  factors  for  the  other  eigenfuntions. 
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As  our  local  displacement- fitting  method  we  match  the  crack  opening  displacement, 
-U/j.  on  the  crack  flank  (0*17)  at  the  five  radial  stations  r  *1,  3,  5,  7,  9.  To  do  this 

0  MM 

we  take  the  five  eigenfunctions  for  the  crack  given  by  (9)  with  n«0,  1,  2,  3,  4 
Then  it  follows  from  (10)  and  the  forms  in  (9)  for  n«0  that  we  must  treat  the  5x5 


system 


K.  ♦  d,r  ♦  d,r2  +  d  r3  »dr4  «  *  .1 —  i 

^  1  2  3  •  2<l-vW  r  "  I  ffmff. 


r  «  1.  3,  5.  7,  9. 

where  ^  is  our  displacement-fitting  estimate  of  K(  and  dn  (n 
participation  factors  for  the  other  eigenfunctions. 


•  1-4)  are  estimated 


**The  eigenfunction  associated  with  n*1/2  is  not  as  given  in  I9>  since  er  g^O  on  8*7T  when  n*1/2  in  (9) 
The  eigenfunction  tor  this  eigenvalue  reouires  individual  analysis  and  eventually  simply  reduces  to  a  equal  to  » 
constant  as  the  only  noniero  stress  in  the  plane.  As  a  result  a  g*0  on  8*0  lor  this  case  so  that  in  effect 
this  eigentuntion  is  also  being  matched  by  our  stress-fitting  procedure. 

***Here  the  actual  eigenfunction  associated  with  n  *  1/2  has  u^  “0  on  8*  17  and  so  contributes  nothing 
to  a  crack  opening  displacement  fit  (just  as  the  corresponding  stress  does  not  affect  a  stress  fit).  Further, 
although  <91  does  give  the  correct  eigenfunctions  for  n«3/2,  5/2.  7/2.  9/2.  none  of  these  terms  add  to  ug 
on  8*1T  Hence  m  effect  here  we  are  fitting  the  first  ten  eigenfunctions. 


/  .*  V  -■  s  V  >  V  "j*  V  ’  * 


s  s 


We  now  apply  our  local  stress  fitting  method  to  Problem  1.  From  Problem  1,  from 
(9),  Table  2,  we  have 


W?  a6\  «Kp 


248  r2  _  rJ 

- r  *  - (190  -  46r  ♦  5r2  ♦  — )]. 

315  189  5 


Hence  system  (12)  becomes 

R,  *  d,  *  d2  *  d3  *  d4  K' 

K,  *  3d1  *  3^d2  ♦  9d3  -  9/3d4  ■  3K. 

K,  ♦  5d,  ♦  SJZd2  *  25d3  ♦  25^5^  ■  5K.  (15) 

Kj  ♦  7d1  ♦  7^7d2  ♦  49d3  -  49/7d4  «  7K, 

Kj  +  9dd  ♦  27d2  ♦  8ld3  *  243d4  ■  9K. 

Solving  (15)  yields  ^  «  0.  d1  »  K.  d2  *  d3  «  d4  ■  0;  the  first  of  these  is  the  estimate 
entered  in  Table  3. 

Similarly  we  apply  our  local  displacement-fitting  method  to  Problem  2.  then  our 
stress-fitting  and  displacement-fitting  methods  in  turn  to  Problems  3  and  4.  The 
resulting  estimates  are  entered  in  Table  3. 

TABLE  3:  STRESS  INTENSITY  VALUES,  K. 


Problem  Basis  of  local  Local  fitting  Actual  exact 

fitting  method  estimate  of  K  value  of  K 


lacement 


>9 


Several  comments  are  in  order  at  this  point.  Clearly,  given  knowledge  of  the  exact 
solution  as  here,  it  is  possible  to  modify  either  of  the  two  local  fitting  methods  so 
that  they  do  work.  In  practice,  however,  one  does  not  know  the  answer  in  an  actual 
application,  so  that  one  does  not  have  the  option  of  fixing  up  a  local  fitting 
approach  but  can  only  see  if  things  "look  right”.  In  this  connection,  consider  the 
plot JTn?  <r^(5=0)  for  Problem  1,  which  as  it  happens  coincides  with  «  7T)12tt1t 

I  2(1  -  V)  for  Problem  2  (Fig.  3).  Certainly  here  things  "look  right".  Indeed  one  could 
even  imagine  obtaining  results  from  a  sequence  of  finite  element  grids  at  the  points 
3.  5,  7.  9,  then  at  2.  3,  4.  5,  6.  7,  8,  9,  and  ultimately  at  3/2,  2.  5/2,  7/2,  4,  9/2,  5, 
11/2,  6,  13/2,  7,  15/2,  8,  17/2,  9,  and  still  concluding  K(  =  0  when,  in  actuality,  K(  * 
K.  And  it  is  possible  to  construct  a  problem  which  has  both  its  stress  fit  and  its 
displacement  fit  do  this.  In  sum,  once  any  local  fitting  procedure  has  been  tested 
and  decided  upon,  it  is  then  possible  to  set  up  an  "application"  on  which  the 
procedure  gives  completely  incorrect  results.  Finally  we  remark  that  while  these 
given  problems  are  not  the  easiest  to  solve  without  hindsight  using  other  methods, 
they  are  by  no  means  impossible.  In  fact,  using  the  superposition  method  of  [3]  in 
conjunction  with  the  path  independent  integrals  of  [4]  we  would  judge  that 
reasonable  results  could  be  obtained  in  return  for  some  computational  effort.  More 
precisely,  for  Problem  1  for  instance,  we  estimate  from  the  participation  of  the 
regular  eigenfunctions  in  Table  2  relative  to  that  in  trial  problems  in  [4],  that  K( 
could  be  found  to  within  10%  using  a  uniform  grid  comprised  of  around  1000 
constant-strain  triangle  elements,  to  within  1%  using  grids  with  up  to  5000  such 
elements  and  extrapolation.  Moreover,  these  analyses  could  be  performed  in  a 
completely  systematic  way  without  drawing  on  a  knowledge  of  the  exact  answer  at 
all. 
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ON  THE  CAPRICIOUS  NATURE  OF  ELASTIC  SINGULARITIES 


AE  Chambers  end  GB  Sinclair 

Department  of  Mechanical  Engineering.  Carnegie-Mellon  University 
Pittsburgh.  PA  15213 

Classical  elastic  analysis  of  re-entrant  corners  or  sharp  notches  in  plates  reveals  the 
presence  of  stress  singularities  at  the  notch  vertices  As  a  result,  one  must  make 
inferences  concerning  the  structural  integrity  of  components  containing  such  features  when 
the  stresses  become  infinite  in  response  to  infinitesimal  loads,  clearly  not  a  simple  task  At 
this  time  the  accepted  methodology,  for  the  case  of  cracks  in  brittle  materials  anyway,  is 
to  take  the  coefficient  of  the  singularity  -  the  stress  intensity  factor  -  as  the  parameter 
governing  failure  This  choice  represents  the  basic  tenet  of  linear  elastic  fracture 
mechanics  today 

The  usual  approach  taken  in  the  elastic  asymptotic  analysis  of  re-entrant  comers  is  to 
examine  what  stress  fields,  including  singular  fields,  ere  possible  for  a  given  local  geometry 
‘Here  we  adopt  an  inverse  approach  of  prescribing  e  field  with  bounded  but  nontrivial 
stresses  at  the  vertex  then  asking  what  problem  can  such  fields  occur  in  In  this  way  we 
construct  a  closed  form  solution  with  finite  stresses  for  e  loaded  plate  in  the  shape  of  a 
pac-man  with  a  specific  angle  between  its  stress-free  flanks  This  configuration  then  has 
a  stress  concentration  factor  but  no  stress  intensity  factor.  By  maintaining  the  same 
loading  on  the  outer  edges  of  the  pac-man  yet  perturbing  the  shape  of  the  flanks  about 
their  original  positions,  problems  which  have  participating  singularities  are  generated  these 
last  do  not  readily  admit  to  the  determination  of  close  form  solutions  but  ere  amenable  to 
numerical  analysis  of  demonstrably  high  resolution  using  appropriate  path  independent 
integrals  and  a  superposition  procedure  Since  an  identical  loading  acts  in  all  problems  and 
only  very  minor  changes  in  geometry  occur,  it  is  to  be  expected  that  the  physics' 
responses  would  be  quite  similar.  Thus  the  stress  intensity  factor's  on-off -on-again  type 
behavior  makes  it  difficult  to  envisage  how  it  could  be  the  controlling  damage  parameter  in 
this  sequence  of  problems  These  results,  end  other  like  them  (eg.  references),  suggest 
the  need  for  a  reconsideration  of  the  role  of  the  stress  intensity  factor  in  fractue 
mechanics 

References 

AE  Chambers  and  GB  Sinclair,  On  obtaining  fracture  toughness  values  from  the  Werstue 
/nt.  J.  Fractur e  30,  R11-R15  (1986) 

GB  Sinclair  and  AE  Chambers.  Strength  size  effects  and  fracture  mechanics  whit  does 
the  physical  evidence  say?  Sngng  Fracture  Mech.  to  appear. 


THE  ELASTIC  ANALYSIS  OF  THREE  RE-ENTRANT  CORNERS 


I 

r. 

■r. 

>. 


SM  86-10 


A  E.  Cnambers  and  G.  B.  Sinclair 


June,  1986 


Department  of  Mechanical  Engineering 
Carnegie  Institute  of  Technology 
Carnegie-Mellon  University 
Pittsburgh,  Pennsylvania  15213 


£ 

S: 

£• 


£ 


THE  ELASTIC  ANALYSIS  OF  THREE  RE-ENTRANT  CORNERS 
INTRODUCTION 

Since  many  engineering  structures  are  designed  with  what  are  known  as  re-entrant 
corners  or  notches  Figure  1),  it  is  important  to  understand  the  mechanical  behavior  of 
these  geometries  upon  loading  Elastic  analysis  of  these  notches  can  result  in  stress  fields 
near  the  notch  tip  that  are  singular  in  nature  1 1  ].  Physically,  these  singularities  do  not 
make  any  sense.  Unbounded  stresses  imply  failure  under  even  an  infinitesimal  load,  which 
is  inconsistent  with  everyday  experience.  For  analytical  purposes,  it  is  desirable  to 
determine  the  participation  of  these  singular  fields  in  comparison  to  the  regular  fields 
present 

These  singular  stress  fields  are  eigenfunctions  of  the  equilibrium  problem  posed  upon 
loading  re-entrant  configurations.  In  general,  most  regular  eigenfunctions  stresses  acting 
along  the  bisector  of  the  re-entrant  corner  are  equal  to  zero,  leaving  only  the  singular 
fields  to  dictate  the  mechanical  behavior  of  the  notch  Figure  11.  It  is  these  stress  fields 
acting  at  the  notch  tip  that  are  analytically  unbounded.  It  is  also  these  stress  fields  that 
require  making  some  sense  of  their  significance  in  the  problem  under  investigation.  In 
order  to  gain  some  understanding  as  to  how  the  singularity  plays  a  role  in  the  overall 
stress  field,  one  could  ask  the  question,  is  there  a  re-entrant  geometry  which  has  non¬ 
zero  and  finite  stresses  acting  at  the  notch  tip?  Such  a  configuration  would  have  no 
singularity  participation,  Le.  a  stress  intensity  factor  of  zero.  The  answer  to  this  question  is, 
in  fact  yes  Upon  finding  a  re-entrant  corner  that  has  only  regular  stress  fields,  the 
objective  of  the  present  work  is  to  track  the  variations  of  the  stress  intensity  factor  with 
alight  perturbations  of  geometry.  A  re-entrant  corner  under  a  particular  loading  that  has  a 
stress  intensity  factor  of  zero  yet  nontrivial  stresses  at  the  vertex  will  be  determined.  The 
angle  near  the  vertex  of  the  corner  will  be  decreased  and  then  increased  slightly.  The 


effects  of  these  geometric  changes  will  be  studied  through  the  behavior  of  the  stress 
intensity  factor. 

The  first  section  presents  the  formulation  of  the  problem  under  investigation.  The 
exact  geometries  to  be  examined  are  introduced.  The  generalized  stress  intensity  factor  is 
defined.  The  second  section  presents  the  numerical  analysis.  A  finite  element  approach  is 
taken  to  determine  the  stress  and  displacement  fields  A  path-independent  integral  is  used 
to  calculate  the  stress  intensity  factor  with  the  information  produced  from  the  finite 
element  analysis.  Results  will  be  presented  and  discussed  in  the  final  sectioa 
PROBLEM  FORMULATION 

In  order  to  study  the  effects  of  a  stress  singularity  of  a  re-entrant  corner,  three 
individual  regions  are  analyzed.  The  basic  difference  of  the  three  regions  is  that  the 
material  angle  local  to  the  notch  tip  is  slightly  perturbed  from  an  initially  symmetric 
geometry.  The  planar  region  R.  has  material  angle,  2a.  (i= 0,1,2).  In  particular,  aQ  - 
128.727°  .  «,  =  125  863  and  «2  =  131.863°  .  The  choice  of  these  particular  values  will 
be  explained  later.  Let  (x,y|  be  cartesian  coordinates  that  have  origin  at  the  notch  vertex, 
and  t^e  x-axis  bisects  the  material  angle  there  (Figure  2).  Let  (r.0)  be  polar  coordinates  in 
the  region  where  r2  =  x2  ♦  y2  and  6  ~  arctan  (y/x).  The  dimensions  of  the  notch  are  as 
follows:  a  =  notch  depth  and  h  *  half-height  of  notch. 

The  region,  RQ,  is  defined  as  the  region  enclosed  by  basically  two  boundaries.  The 
first  boundary.  8 1RQ,  is  the  boundary  to  which  the  prescribed  tractions  will  be  applied  It 
consists  of  three  sides  of  a  rectangle.  The  line  segments  that  define  it  are. 

y  =  ±h  for  -a£x£  2h  and  x  *  2h  for  -hSySh  (1) 

The  second  boundary.  92R0-  to  *001086  the  region,  RQ,  is  the  two  line  segments  given  by 

tan«0  *  ±y/x  for  -a£x£0.  (2) 

The  region,  R^  is  defined  in  a  similar  manner.  The  local  (x,y)  ci-O’dinate  system  near 
the  vertex  of  the  corner  is  slightly  rotated  about  the  global  coordinates  by  about  5  ~  3°. 
The  outer  rectangular  boundary,  8^  .  on  which  the  tractions  are  applied  is  left  intact,  as 


is  the  lower  flank.  The  upper  flank  is  rotated  by  approximately  23  and  extended  to  y  = 
h/2,  thus  resulting  in  less  solid  material  to  support  the  applied  loads  in  the  original  problem. 
It  is  connected  by  a  smooth  curve  to  the  point  where  the  positive  y-face  boundary  begins 
in  the  second  quadrant  at  y  =  h  This  is  boundary  32^i 

The  perturbation  of  the  upper  flank  for  Problem  i  =  2  is  similar  to  that  in  the  first 
perturbed  problem  However,  the  rotation  of  the  local  axes  is  in  the  opposite  direction, 
resulting  in  more  solid  material  than  both  Problems  i  =  0  and  1.  This  is  region  R  (The 
dotted  lines  in  Figure  2  represent  the  upper  flanks  for  Problems  i  =  1  and  2.) 

Formally,  we  seek  the  stresses,  arf  .  a  gg  and  rfg  .  and  the  displacements,  Ug  and  uf, 
as  funcjions  of  r  and  8  in  the  region  R.,  satisfying  the  following.  The  stress  field  must 
satisfy  the  two-  dimensional  stress  equations  of  equilibrium  in  the  region  in  the  absence  of 
body  forces. 

da  1  dr  a  a  ~aaa 

_ u  + _ -is.  +  .  nr.— 1 88.  =  n 

dr  r  dd  r 


1  A  n  dr  n  2 7  n 

_ vv  +  .  r&  +  -  r.B  =  o, 

r  d8  dr  r 


(3) 


on  R(.  The  stress-displacement  relations  for  a  homogeneous,  isotropic,  linear  elastic  plane 
hold  for  the  regions. 
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on  R.,  where  //  is  the  shear  modulus,  C1  equals  1-2v  for  plane  strain  and  1-v  for  plane 
stress,  C2  equals  1-v  for  plane  strain  and  1  for  plane  stress,  and  v  is  Poisson's  ratio 
The  notch  flanks  are  stress-free  which  requires  that 

o qq  =  r«  =  0  on  B2  R.  when  6  =  ±  a.  (5) 

The  remaining  boundary  conditions  on  8^  are  prescribed  tractions  derived  from  the 
following  stresses: 
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where  <?0  is  arbitrary.  Finally,  the  regularity  requirements  involve  the  fact  that  the 
displacements  at  the  vertex  of  the  notch  must  be  finite, 

u  and  ua  -  0(1)  on  R  as  r-»0.  (7) 

r  o  i 


More  specifically,  we  want  to  use  the  stress  and  displacement  fields  found  in  the 
problems  (i  =  0,  1  and  2)  formulated  above  to  find  the  generalized  stress  intensity  factors 
present  There  are  two  possible  factors,  one  arising  from  any  symmetric  loading  and  one 
arising  from  any  antisymmetric  loading  in  the  problem.  For  the  symmetric  case,  the  stress 
intensity  factor  is  defined  as: 

K,  =  oeQ  r1'*.  on  0  =  0.  (8) 

The  subscript  I  represents  mode  I  loading  (symmetric  opening  mode).  The  value  of  X  is 
the  eigenvalue  identified  after  Williams  [  1  ]  that  represents  the  only  singular  stress  field 
possible  at  the  re-entrant  corner  under  this  symmetric  mode.  That  is,  X.  is  the  root  of  the 
eigen-condition 

sin2X  a  =  -X  sin2«  (0<X  <1).  (9) 

ii  ii  i 


For  the  anti-symmetric  case,  the  stress  intensity  factor  is  defined  as: 
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where  the  subscript  II  represents  the  mode  II  loading  (antisymmetric), 
evaluated  from  the  antisymmetric  eigen-condition, 


The  eigenvalue 


(10) 

X  is 
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K 
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sin2A.<r,  =  A.  sin2«.  (0<A.<1).  (11) 

These  eigenvalues  produce  eigenfunctions  which  are  the  only  singular  stress  fields  possible 
for  the  particular  notch  geometry.  As  a  result  of  the  problem  formulation,  the 
displacement  fields  derived  are  regular.  That  is,  the  fields  near  the  notch  tip  behave  as 

oqq  and  r =  0(r\  1  )  and  uf  and  u^  =  0(r\  )  on  R>  as  R-»0.  (12) 

Some  explanation  of  the  traction  boundary  conditions  found  from  (6)  is  now  in  order. 
Recall,  that  the  purpose  of  this  study  is  to  look  at  slight  geometric  perturbations  of  a 
regular  configuration.  Therefore,  RQ  is  taken  to  be  the  region  in  the  nonsingular  problem 
under  investigation.  The  approach  taken  is  such  that  a  particular  material  angle,  aQ,  is 
sought  that  will  produce  a  regular  and  non-zero  stress  field  for  an  antisymmetric  problem. 
The  following  transcendental  equation  results. 

tan2«0  =  2«0  (13) 

It  is  in  solving  this  equation  that  aQ  is  found  to  be  approximately  128.7°  .  The  remaining 
two  problems  are  then  posed  by  simply  changing  the  material  angle  by  plus  or  minus  about 
three  degrees.  Note,  also,  there  is  now  an  exact  solution  for  the  nonsingular  problem  that 
is  independent  of  the  radial  coordinate,  r.  It  is  these  regular  stress  fields  that  are  used  as 
the  traction  conditions  on  the  boundary  8^.  (i=0,1,2). 

In  conclusion,  what  we  have  are  three  problems  where  the  applied  tractions  and  the 
boundary  on  which  they  are  applied  remain  the  same  The  flanks  change  in  geometry,  but 
remain  stress  free  An  easy  way  to  imagine  it  all,  is  that  the  region  Rp  is  an  ideal 
geometry  of  a  re-entrant  corner  that  is  to  be  fabricated  (i=0).  The  first  singular  case 
examined  (i-1)  is  a  result  of  machining  just  a  little  too  much  material  away  from  the  upper 
flank  The  second  (i=2)  follows  from  not  machining  quite  enough  material  away  from  the 
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flank.  When  the  three  configurations  are  loaded  under  the  same  conditions,  what  is  the 
result  of  these  slight  geometric  changes? 
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ANALYSIS 

The  analysis  of  the  problems  requires  two  major  steps.  First,  the  stress  fields  in  R  (i 
=  0,1,2  )  need  to  be  determined.  Secondly,  as  stated  in  the  problem  formulation  section,  a 
method  is  then  needed  to  determine  the  stress  intensity  factor  of  each  problem  for  both 
mode  I  and  mode  II  loadings. 

A  finite  element  method  is  used  to  numerically  determine  the  resultant  stress  and 
displacement  fields  of  the  three  problems*  For  each  problem,  three  grids  of  different 
sizes  are  used  in  the  hopes  of  attaining  better  accuracy  and  gauging  convergence.  The 
grid  refinement  is  fairly  systematic.  The  coarse  grid  for  the  nonsingular  problem  (  i=0  )  is 
shown  in  Figure  3  The  nonsingular  stress  fields  defined  in  (6)  are  applied  on  the 
rectangular  boundary,  leaving  the  notch  flanks  stress-free.  Upon  finite  element  analysis, 
nodes  along  common  8  are  of  constant  stress  (One  drawback  of  the  grid  pattern  chosen 
is  that  near  the  notch  vertex  there  is  no  refinement  in  the  angular  direction.  Therefore, 
convergence  as  r  goes  to  zero  is  difficult  to  judge,  though,  divergence  can  be  obvious) 
The  coarse  grid  for  the  first  singular  problem  (i=1)  is  shown  in  Figure  4,  while  that  for  the 
second  singular  problem  is  displayed  in  Figure  5.  The  dotted  lines  in  Figures  4  and  5  are 
the  smooth  curves  that  represent  the  flanks.  Since  the  finite  element  code  uses  linear 
elements,  these  curves  are  better  approximated  upon  grid  refinement  For  all  three 
problems,  three  grids  were  used  in  the  stress  analysis.  The  coarse  grids  have  68  elements 
and  48  nodes.  The  medium  grids  have  264  elements  and  159  nodes  The  fine  grids  have 
1040  elements  and  573  nodes. 

The  second  step  of  the  analysis  requires  a  method  of  determining  the  generalized 
stress  intensity  factor  from  the  stress  and  displacement  fields  generated  from  the  finite 
element  analysis  The  method  employed,  the  H-integral  [2],  is  from  the  family  of  path- 
independent  integrals.  The  major  advantage  of  path-independent  integrals  is  that  the  only 

The  program  used  is  PlANDJ,  originally  developed  by  J.L.  Swedlow.  It  is  used  on  the  VMS/VAX  at  the 
Mechanical  Engineering  Department.  Carnegie-Mellon  University. 
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errors  accumulated  are  those  from  the  numerical  approximations  present  Another 
advantage  of  these  integrals  is  the  fact  that  the  symmetric  and  antisymmetric  parts  need 
not  be  separated  in  the  analysis.  (For  the  argument  supporting  this  statement,  see 
reference  [33.)  The  path  of  integration  is  shown  in  Figures  3,  4  and  5.  Notice  that  the 
inside  path  is  chosen  so  that  the  curvature  of  the  notch  flanks  in  the  two  singular 
problems  (i  =  1  and  2)  does  not  interfere  with  the  independence  of  the  integration  path. 

The  H-integral  program  uses  the  second-order  trapezoidal  rule  for  the  integration.  The 

nodal  values  along  the  path  of  integration  are  used  as  points  of  numerical  integration.  The 
complementary  stress  fields  and  displacement  fields  used  in  the  H-integral  (see  reference 
[2])  are  determined  from  the  (r .6)  coordinates  of  the  nodes  along  the  path.  The  actual 
nodal  displacements  are  read  from  the  finite  element  output  The  stresses  at  the  nodes  are 

determined  from  nodal  averages  of  the  surrounding  elements  of  each  node  along  the 

integration  path. 

RESULTS 

In  presenting  the  results  of  this  analysis,  some  sort  of  verification  is  required.  One 
common  method  of  verification  of  the  finite  element  method  is  to  check  convergence  at 
nodes  common  to  all  three  grids  used  for  each  problem.  In  the  nonsingular  case, 
convergence  to  the  exact  solution  (6)  is  required.  In  the  other  two  cases,  convergence 
within  the  three  associated  grids  is  desired.  As  stated  earlier,  since  there  is  no  grid 

refinement  in  the  angular  direction,  convergence  near  the  origin  is  difficult  to  assess. 

Divergence,  on  the  other  hand,  is  clear  when  present  Therefore,  the  nodal  stresses  acting 

along  the  x-axis  (0=0)  are  plotted  for  alt  three  problems.  (Figure  6,  7,  and  8.) 

Convergence  to  the  exact  solution  in  the  nonsingular  problem  (i=0)  is  seen  clearly  (Figure  6). 
The  antisymmetry  in  the  problem  is  observed  in  that  a qq  ~  0  and  r ^  is  finite  along  the 
x-axis  Convergence  at  r=0  for  r r0  is  not  obvious  (though,  as  explained  earlier,  it  is  not 
particularly  expected).  However,  divergence  is  not  noticed. 

Comparing  the  expected  singular  problems  (i  =  1  and  2).  divergence  at  the  origin  is 
clear.  For  the  singular  problem  i=1  Figure  7),  a qq  at  r  =  0  is  divergent,  where  the  trend 
of  T'g  is  not  obvious  This  is  consistent  with  what  will  be  seen  later,  problem  i  =  1  has  a 


R avrc,  (o  • 


mode  I  singularity  and  no  mode  II.  The  results  for  and  r fg  for  problem  i=2  (Figure  8) 
are  divergent  at  r=0.  This  is  also  consistent  since  both  mode  I  and  II  singularities  are 
present  in  the  problem. 


Another  method  of  verification  is  the  fact  that  where  there  is  no  singularity  present, 
H  =  K  should  be  zero  The  H-integral  of  the  nonsingular  problem  oscillates  about  zero. 
And  is  in  fact,  numerically  zero  for  our  purposes.  (The  H^-integral  cannot  be  determined 
for  the  e'genva'ue  equal  to  1.0,  see  reference  [2].) 

For  the  nonsingular  problem  (i=0),  the  eigenvalues  corresponding  to  both  equations  (9) 
and  (11)  are  simply  equal  to  1.0.  Therefore,  the  stress  and  displacement  fields  are  of  the 
order  one  The  mode  I  eigenvalue  of  the  first  singular  problem  (i=1)  is  X  1  =  0.5818  this 
is  the  only  singular  stress  field  in  this  problem.  For  the  second  singular  problem  (i=2). 
however,  singularities  in  both  modes  exist  That  is,  X^  =  0.5554  and  X2n  =  0.9528.  The 
corresponding  results  of  the  H-integration  are  presented  in  Table  I. 

Notice  how  the  value  of  the  dimensionless  in  problems  i=1  and  2  are  virtually  equal 
in  magnitude  but  opposite  in  sign.  This  could  have  been  expected  from  the  beginning  since 
the  original  (nonsingular)  material  angle  is  perturbed  about  the  same  in  both  directions  for 
each  of  the  singular  problems  (i.e.  a1  =  aQ-$  and  «2  =  «0+$).  Notice,  also,  that  the  value 
of  K(|  for  problem  i=2  is  larger  in  magnitude  than  K(.  This  is  expected  since  the  loading 
used  was  derived  from  an  antisymmetric  problem.  The  antisymmetric  singularity  should  then 
be  dominant 

CONCLUDING  REMARKS 

In  general,  re-entrant  corners  produce  stress  fields  which  are  singular  as  the  notch 
vertex  is  approached.  However,  a  configuration  can  be  found  at  a  particular  material  angle 
and  loading  that  will  produce  all  regular  stress  fields  Slight  perturbations  in  geometry 
produce  predictable  behavior  in  the  stress  intensity  factor.  Since  the  applied  loading  in  this 
particular  problem  is  based  on  an  antisymmetric  configuration,  the  mode  II  singularity  has  a 
larger  participation  factor  than  that  of  the  mode  I.  Perturbing  the  material  angle  in  equal 
magnitude  but  opposite  direction,  results  in  equal  but  opposite  stress  intensity  factors  for 
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Table  I.  Dimensionless  Stress  Intensity  Factors  (K  *  K  /  #1  M 


Problem/ 

Mode 


Problem  0 


Coarse  Medium  Fine 


|  Mode  1 

1.000 

128.7° 

0000 

0.000 

0.000 

Problem  1 

Mode  1 

0582 

1259° 

-0.066 

-0060 

-0  058 

Problem  2 

}  Mode  1 

0555 

131.9“ 

0.063 

0.060 

0059 

f  Mode  II 

0.953 

131.9° 

0.648 

0.610 

0632 

X  =  eigenvalue  after  Williams  [1] 
a  -  half  the  material  angle  of  re-entrant  corner 


Even  though  the  behavior  of  the  stress  intensity  factor  may  be  predictable,  the 
problem  still  exists  in  that  these  singular  fields  are  not  able  to  interact  with  the  other 
regular  stress  fields  that  are  present  The  divergent  stresses  acting  along  y  =  0  in  the 
singular  problems  (i  =  1  and  2)  highlight  this  difficulty.  Studying  the  behavior  of  the  stress 
intensity  factor  may  lead  to  a  better  understanding  of  what  is  occuring  in  these  re-entrant 
configurations 
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ON  OBTAINING  FRACTURE  TOUGHNESS  VALUES  FROM  THE  LITERATURE 
A.E.  Chambers  and  G.B.  Sinclair 
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Recently,  extensive  compendlums  of  sources  of  fracture  toughness  data 
have  been  assembled  by  Hudson  and  Seward  [1,2].  These  lists  of  references 
naturally  prompt  the  question  as  to  what  an  engineer  would  find  upon  con¬ 
sulting  them  in  terms  of  actual  values  of  fracture  toughness  for  a  given 
material.  The  intent  of  this  note  is  to  indicate  an  answer. 

In  choosing  specific  metals  for  which  to  obtain  fracture  toughness 
values,  we  look  to  ones  with  large  listings  of  data  sources  so  as  to 
gauge  any  variability  present.  To  this  end  we  select  AISI  4340  steel  and 
7073-T6  aluminum.  We  focus  on  plane  strain  fracture  toughness  values 
(K^)  as  governed  by  ASTM  standards  since  these  are  more  generally  regard¬ 
ed  as  being  material  properties.  We  do  not  check,  however,  whether  or  not 
a  given  test  furnishing  a  K  value  complies  with  all  of  the  specifications 
in  ASTM  E399  [3],  simply  because  none  of  the  references  reviewed  provided 
sufficient  information  to  enable  a  complete  check.  Hence  values  are  in¬ 
cluded  as  valid  K_  if  their  contributors  claim  them  as  such.  Since  we 
frequently  encountered  plane  stress  fracture  toughness  values  (K  )  in  the 
data  for  7075-T6  aluminum,  we  include  these  as  a  separate  set  for  com¬ 
parison.  Further  by  way  of  comparison,  we  note  yield  strengths  (o  ), 
since  a  may  reasonably  be  viewed  as  the  uniaxial  tension  test  quantity 
analogous  to  initial  unstable  crack  propagation  in  a  brittle  material. 

In  processing  the  data  we  distinguish  between  markedly  different  specimen 
types  reported  in  a  single  source  but  otherwise  use  mean  values  for  each 
source.  That  is,  when  what  is  in  essence  the  same  test  is  repeated  a 
number  of  times  and  outcomes  recorded,  we  merely  extract  the  mean*.  For 
these  average  values,  we  note  the  number  of  essentially  Independent  sources 
and  calculate  an  overall  mean.  To  estimate  the  variability  we  also  cal¬ 
culate  the  ranges  and  95  percent  confidence  intervals  (1.96s,  8  being  the 
standard  deviation).  In  justification  of  the'  second,  histograms  of  the 
data  show  good  agreement  with  the  expected  frequencies  in  a  normal  distri¬ 
bution  (possibly  because  the  individual  data  typically  represent  mean 
values  themselves  and  the  central  limit  theorem  applies  to  some  extent). 

The  only  exception  to  this  agreement  occurs  for  the  K_  values  for  7075-T6. 
These  data,  though,  conformed  well  with  a  normal  distribution  on  taking 
logs  and  accordingly  a  log  transformation  was  used  to  determine  the  confi¬ 
dence  Interval  in  this  instance.  The  results  are  summarized  in  Table  1; 
details  of  K  values  are  shown  in  Tables  2,3  wherein  single  numbers  in 
brackets  denote  original  sources,  hyphenated  numbers  the  corresponding 
reference  in  [1]  or  [2],  and  a  virgule  between  the  two  implies  as  reported 
in  the  latter,  e.g.  [6]/[l-12]  is  [6]'s  data  as  drawn  from  [12]  in  [1]**. 
*For  the  7075-T6  K  ,  about  a  quarter  of  the  results  were  designated  as 
being  for  a  longitudinal  orientation  and  a  like  fraction  as  being  for  a 
transverse  orientation  with  the  remainder  being  unspecified.  Given  that 
the  difference  between  the  means  for  all  the  longitudinal  and  transverse 
cases  was  found  to  be  less  than  2.5  percent,  the  effects  of  orientation 
for  this  alloy  were  not  regarded  as  being  sufficiently  significant  to  merit 
distinct  classification. 

**In  the  interests  of  brevity  we  do  not  relist  references  given  in  [1],  [2]. 
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Examining  Table  1,  we  find  that  the  variability  in  plane  strain  frac¬ 
ture  toughness  for  AISI  4340  steel  is  100  percent  of  the  overall  mean  when 
based  on  the  range,  87  percent  when  based  on  the  confidence  interval  The 
corresponding  percentages  for  7075-T6  aluminum  are  162  percent  and  108 
percent.  These  wide  fluctuations  are  in  marked  contrast  to  the  respective 
percentages  for  the  yield  strengths  which  are  17  percent  and  19  percent 
for  the  steel,  and  22  percent  and  21  percent  for  the  aluminum.  Indeed 
they  are  quite  comparable  to  the  variations  in  plane  stress  fracture 
toughness  values  for  7075-T6,  viz,  110  percent  and  72  percent*. 

Such  diverse  K.  data  demonstrate  that  the  engineer  needs  to  exercise 
considerable  care  in  drawing  a  value  of  plane  strain  fracture  toughness 
for  a  particular  material  from  the  literature.  Short  of  undertaking  the 
time  consuming  task  of  collecting  a  sufficient  set  of  references  as  a  way 
of  assessing  variations  as  here,  no  reasonably  certain  methodology  for 
ensuring  a  conservative  estimate  of  fracture  toughness  appears  to  exist. 
Possibly,  if  a  single  value  is  found  and  then  assumed  to  be  at  the  upper 
limit  of  the  ranges  involved  and  a  safety  factor  applied  to  reduce  it  to 
the  lower,  a  conservative  could  be  expected.  Such  safety  factors 
here  would  be  2.9,  2.5  for  ASsi  4340  and  3.7,  2.9  for  7075-T6.  These 
suggest  that  a  safety  factor  of  3  might  be  adequate,  although  there  is 
really  no  guarantee  of  this  being  so  for  another  material.  Indeed,  the 
fact  that  the  plane  strain  fracture  toughness  data  display  more  than  five 
times  greater  variability  than  the  yield  stress  data  and  vary  about  as 
much  as  the  plane  stress  values  which  are  known  to  be  geometry  dependent, 
raises  serious  doubts  concerning  the  notion  of  KIc  being  a  material  pro¬ 
perty.  An  engineer  continuing  to  interpret  and  use  it  as  such  may  there¬ 
fore  be  staking  significant  errors. 
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Material  Quantity  No.  of  Mean 

data  value 


AISI  4340 


Range 


5% 

Confidence 

interval 


1370-1630 

1360-1650 

37.4-108.2 

39.8-101.2 

434-545 

452-561 

21.0-76.9 

19.4-56.7 

32.9-111.0 

45.1-96.2 

Aluminum- 
707  5-T6 


Note:  O'  values  are  in  MPa;  K_  ,  K  values  are  in  MPa  m 
y  Ic  c 

Table  2.  Plane  strain  fracture  toughness  values  for  AISI4340  steel 


K 

2  Ic 
)  Source  (Mpa  m 

11 2 )  Source 

[111/11-15] 

[1-25] 


[1-152] 


[12]/ [1-156] 
[1-160] 

[  13 ] / [ 2—3 ] 


[14) /[2-3] 
[2-27] 

[15] / [2-130] 
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Table  3.  Plane  strain  fracture  toughness  values  for  7075-T6  aluminum 


(MPa  m1/2) 

Source 

Kt 

lc  1/2 
(MPa  m  '  ) 

Source 

(MPa  m1^2* 

Source 

37.0 

[161/1-2] 

31.6 

[20] / [1-51] 

26.6 

[2-66] 

21.0 

[1-5] 

30.8 

[1-52] 

28.2 

[24] / [2-130] 

27.7 

24.2 

[2-2] 

31.9 

[25] / [2-130] 

38.5 

[1-6] 

39.1 

[ 21 ] / [ 2— 3 ] 

29.9 

[26] / [2-130] 

48.9 

25.3 

[22] / [2-3] 

26.4 

[27] / [2-137] 

35.9 

[1-11] 

34.5 

[23]/2-27] 

35.2 

[2-153] 

48.1 

76.9 

21.6 

[2-157] 

53,7 

34.8 

[2-39] 

32.6 

[28]/ [2-159] 

29.1 

[17]/ [1-12] 

35.7 

[2-59] 

34.8 

[29] / [2-159] 

28.6 

[18 ] / ] 1— 15 ] 

30.8 

[19]/ [1-24] 
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